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Executive  Summary 


Objectives: 

A  new  multiple  access  method,  called  Collision  Frequency  Shift  Keying  (CFSK),  is  investigated  for 
wireless  communication  systems.  The  specific  objectives  of  this  work  are  to: 

•  Develop  a  reduced  complexity  iterative  receiver  for  CFSK  systems. 

•  Simulate  the  performance  of  highly  loaded  CFSK  communication  systems  using  iterative 
receivers. 

•  Perform  thorough  performance  comparisons  of  CFSK  and  CDMA  systems. 

•  Design  and  test  detectors  suitable  for  analog  implementations  of  iterative  CFSK  receivers. 

Accomplishm  ents; 

The  theoretical  capacity  of  the  multiuser  CFSK  system  is  well  known  to  come  at  the  expense  of  an 
optimal  multiuser  receiver  whose  complexity  is  exponential  in  the  number  of  the  users.  In  fact,  we 
conjecture  that  the  optimal  detection  problem  is  provably  NP  complete.  In  order  to  achieve  the 
theoretical  promise  of  the  CFSK  system  and  to  improve  on  the  capacity  of  a  comparable  CDMA  system  a 
reduced  complexity,  suboptimal  receiver  is  necessary.  We  have  developed  two  novel  iterative  receivers. 

The  iterative  receivers  are  modifications  of  the  maximum  a  posteriori  (MAP)  based  consensus  decoder. 
The  metric  of  the  MAP  algorithm  in  the  noisy  CFSK  channel  is,  in  general,  too  computationally  complex 
in  a  highly  loaded  environment.  To  reduce  this  complexity  to  a  usable  level,  we  have  developed  the 
simplified  K-Most  Probable  Combination  (KMPC)  and  Two-Most  Probable  Combination  (TMP)  metrics. 
Both  of  these  metrics  substantially  reduce  the  complexity  of  the  overall  receiver.  Simulation  results 
have  shown  that  this  metric  results  in  a  modest  loss  of  performance  compared  to  the  optimal  receiver. 

Simulation  results  show  that  the  CFSK  system  with  our  iterative  receivers  outperform  conventional 
CDMA  systems.  These  results  strongly  indicate  that  CFSK  is  a  promising  system  for  ad  hoc 
communications  networks  where  a  highly  coordinated  system  like  CDMA  is  not  practical.  The  iterative 
receivers  utilize  MAP  decoders  that  have  been  showm  to  have  viable  low  power,  high  speed  analog 
implementations. 
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Chapter  1 
Introduction 

In  the  last  several  decades,  researchers  have  witnessed  rapid  developments 
and  widespread  applications  in  the  area  of  multiple  access  communication 
systems,  especially  of  mobile  communication  systems.  In  the  first  genera¬ 
tion  of  mobile  communication  systems,  analog  transmission  using  frequency 
division  multiple  access  (FDMA)  [1]  was  used.  In  such  systems,  each  user 
occupies  a  separate  frequency  band  for  transmission.  When  a  large  num¬ 
ber  of  users  access  the  system,  a  huge  amount  of  frequency  resources  are 
required.  In  other  wireless  multiple  access  communication  systems,  such  as 
the  European  recommendation  for  Globe  Speciale  Mobile  (GSM)  communi¬ 
cation  system  [2],  time  division  multiple  access  (TDMA)  is  used.  For  TDMA, 
all  users  share  the  same  frequency  band  by  dividing  the  channel  into  small 
time  slots  and  assigning  the  time  slots  to  different  users.  However,  both 
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FDMA  and  TDMA  are  not  efficient  in  utilizing  the  channel  resoiurces.  Huge 
bandwidth  is  required  for  both  systems  when  a  lot  of  users  access  them. 

Spread  spectrum  communication  originated  in  military  appHcations  in  the 
1950s.  Through  the  appearance  of  the  famous  IS95  standard  proposed  by 
Qualcomm  [3]  [4]  [5],  more  people  realized  that  spread  spectrum  multiple  ac¬ 
cess  (SSMA)  is  much  more  efficient  than  TDMA  and  FDMA.  SSMA  is  also 
known  as  code  division  multiple  access  (CDMA)  [6],  where  pseudonoise  (PN) 
sequences  are  used  to  transform  the  signals  from  other  users  into  a  Gaussian 
like  noise  signal.  The  PN  sequence  can  be  applied  in  the  time  domain,  which 
is  called  direct-sequence  SSMA  (DS-SSMA),  or  in  the  frequency  domain, 
which  is  called  frequency-hop  SSMA  (FH-SSMA).  Each  user  in  these  sys¬ 
tems  is  uniquely  identified  by  their  PN  code.  Thus,  the  signals  are  hopping 
in  the  time  or  firequency  domain  with  a  hopping  pattern  controlled  by  the  PN 
code.  The  conventional  detection  method  for  the  DS-  and  FH-SSMA  systems 
is  not  optimal  in  the  sense  that  it  implements  single  user  detection  by  treat¬ 
ing  other  users  as  purely  unwanted  interference.  Optimal  detectors  require 
multiuser  detection  techniques  that  treat  other  users’  information  as  useful 
in  detecting  a  user’s  information  sequence.  A  variety  of  multiuser  detec¬ 
tion  methods  for  DS-SSMA  systems  [7]  [8]  [9]  [10]  and  for  FH-SSMA  systems 
[11]  [12]  have  been  proposed  in  resent  years. 
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The  third  category  of  multiple  access  systems  is  based  on  collision-type  chan¬ 
nels.  Network  communications  based  on  the  ALOHA  protocol  [13]  is  a  fa¬ 
mous  example  of  the  collision-type  channel.  Communications  on  collision 
type  channels  is  based  on  the  avoidance  of  collisions.  It  is  not  difficult  to 
show  that  the  ALOHA  system  is  even  less  efficient  than  the  TDMA  system. 
However,  resent  research  on  trellis-coded  multiple  access  (TCMA)  gives  us  a 
new  approach  for  communication  on  collision-type  channels  [14].  For  TCMA, 
all  users  use  trelhs-coded  modulation  schemes  and  share  a  common  signal 
space  for  the  modulation.  Thus,  when  different  users  transmit  on  the  same 
signal  a  collision  occurs.  In  [15],  Brannstrdm  showed  that  the  TCMA  channel 
is  nior6  0ffici6iit  tlisiii  tli6  CDM[A  systBiii. 

This  research  project  considered  a  new  multiple  access  communication  system 
with  multiuser  detection  called  the  collision  frequency  shift-keyed  (CFSK) 
system.  The  CFSK  channel  is  a  kind  of  collision  type  multiple  access  chan¬ 
nel  [16]  where  different  users  share  the  common  channel  simultaneously  and 
thus  collisions  occur.  For  each  user,  M-ary  frequency  shift-keyed  (FSK)  [6] 
modulation  is  used  and,  since  all  users  share  the  same  set  of  frequencies, 
collisions  happen  when  different  users  transmit  the  same  frequency  at  the 
same  time.  Hence,  this  is  different  from  an  FDMA  system  where  each  user 
IS  assigned  a  subband  of  the  whole  spectrum  and  no.  collisions  happen  on 
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the  separated  subbands.  In  a  CFSK  system,  users  send  their  signals  within 
the  same  wideband  spectnun,  which  is  similar  to  the  CDMA  system  where 
each  user’s  channel  signal  will  span  the  whole  wideband  spectrum.  How¬ 
ever,  in  the  CFSK  system  one  user  only  transmits  a  single  frequency  point 
in  a  symbol  period.  In  a  CDMA  system,  a  signal  always  occupies  the  whole 
spectrum. 

There  are  some  similarities  and  differences  between  the  CFSK  system  and 
a  FH-SSMA  system.  In  both  systems,  the  users  may  transmit  the  same 
frequency  simultaneously  and  a  collision  happens.  Also,  the  transmitted 
frequencies  in  different  symbol  periods  are  changing  in  a  random-like  pattern. 
In  the  CFSK  system,  the  frequency  hopping  pattern  is  controlled  by  each 
user’s  channel  code,  and  thus  controlled  by  each  user’s  information  sequence. 
On  the  other  hand,  in  an  FH-SSMA  system  the  frequency  hopping  code  is 
given  by  a  PN  sequence  that  is  not  dependent  on  the  user’s  information 
sequence. 

Chang  and  Wolf  first  derived  the  channel  capacity  for  the  noiseless  T-user 
M-frequency  multiple-access  channel  with  and  without  intensity  information 
[17]  in  1982.  In  [18],  Grant  and  Schlegel  compared  the  results  of  Chang  and 
Wolf  and  the  capacity  of  the  CDMA  system  given  by  [19]  [20]  [21].  The  results 
show  that  the  capacity  for  the  CFSK  channel  with  intensity  information  has 
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a  higher  sum  channel  capacity  than  the  conventional  and  multiuser  CDMA 
systems,  especially  when  the  spreading  factor  is  of  moderate  values.  In  Figure 
1.1,  the  sum  capacity  for  the  synchronous  CPSK  system  and  CDMA  system 
are  shown.  When  the  spreading  factor  is  of  moderate  or  high  values,  the 
sum  capacities  of  the  CFSK  channel  with  or  without  intensity  information 
are  almost  4  times  the  sum  capacity  of  the  CDMA  system  with  single  user  or 
multiuser  detection.  When  the  spreading  factor  is  lower  than  one,  the  CFSK 
channel  with  intensity  information  still  has  a  higher  sum  capacity  than  that 
of  the  CDMA  system  with  multiuser  detection.  However,  the  sum  capacity 
for  the  CFSK  channel  without  intensity  information  decreases  drastically 
when  the  spreading  factor  is  lower  than  1.  The  study  of  the  CFSK  system 
with  intensity  information  when  the  spreading  factor  is  of  moderate  values 
is  of  considerable  interest.  This  is  the  emphasis  of  the  work  in  this  project. 

For  asynchronous  CFSK  channels,  the  spectral  efficiency  decreases,  especially 
when  the  number  of  users  is  larger  than  the  number  of  frequencies,  i.e.,  the 
spreading  factor  is  lower  than  1.  The  performance  of  the  asynchronous  CFSK 
system  is  discussed  in  detail  in  Chapter  3. 

There  are  other  advantages  of  the  CFSK  system  over  conventional  multi¬ 
ple  access  systems  which  provide  additional  motivation  for  investigating  the 
CFSK  system.  First,  it  is  very  resistant  to  the  near-far  effect,  as  the  re- 
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Spreading  Factor  =  N/K  (chips  per  user) 


Figure  1.1:  A  comparison  of  the  sum  capacity  of  CDMA  and  CFSK  systems. 
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ceiver  utilizes  multiuser  detection.  Verdu  proved  that  multiuser  detection 
can  provide  good  resistance  to  near-far  effect  in  [7].  Second,  inter-symbol 
interference  (ISI)  is  no  longer  a  problem  in  the  CFSK  system  since  aliasing 
in  the  frequency  domain  is  of  more  significance  than  aliasing  in  time  do¬ 
main.  Third,  because  of  the  frequency  diversity  of  the  CFSK  system,  its 
performance  on  frequency  selective  fading  channels  is  better  than  that  of 
conventional  multiple  access  system.  Finally,  the  transmitter  of  the  CFSK 
system  is  of  very  low  complexity. 

The  channel  metric  for  the  synchronous  and  asynchronous  CFSK  channel 
with  intensity  information  is  given  by  Arnold,  Schlegel  and  Perez  [22]  in 
1996.  Based  on  the  channel  metric,  an  optimal  receiver  can  be  constructed 
using  the  maximum-likelihood  sequence  detection  (MLSD)  method  [7].  Aulin 
and  Espineira  discussed  MLSD  for  the  TCMA  system  in  [14].  The  CFSK 
system  can  in  fact  be  viewed  as  a  multi-dimensional  version  of  a  TCMA 
system  where  each  frequency  spans  a  dimension  of  the  signal  space.  The 
MLSD  method  constructs  a  super-code  and  a  super-trellis  by  considering  all 
users’  memories  and  inputs  together,  and  does  maximum-likelihood  detection 
on  the  super-trellis.  When  the  number  of  users  gets  large,  the  states  of 
the  super-trellis  increase  exponentially  and  the  MLSD  method  is  impossible. 
Thus,  a  suboptimal  method  using  iterative  detection  is  preferred.  This  is 
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another  topic  addressed  in  detail  in  this  project. 

This  report  begins  with  a  description  of  the  CFSK  system  in  Chapter  2. 
The  spectral  efficiency  of  the  synchronous  and  asynchronous  CFSK  channels 
is  discussed  in  Chapter  3,  which  provides  a  theoretical  foundation  for  later 
discussions.  Several  iterative  multiuser  detection  techniques  are  discussed 
in  Chapter  4.  The  performance  of  the  iterative  dectectors  is  discussed  in 
Chapter  5.  These  results  show  that  decent  performance  can  be  achieved 
for  the  system  with  low  and  moderate  loads.  In  Chapter  6,  the  issue  of 
code  design  for  the  CFSK  system  is  considered.  Finally,  conclusions  and 
recommendations  for  further  research  are  given  in  Chapter  7. 


Chapter  2 

CFSK  System  Description 

2.1  Overview 

A  multiple-access  multiuser  system  based  on  the  Collision  Frequency  Shift- 
keyed  (CFSK)  channel  is  described  in  Figure  2.1.  Suppose  that  there  are  T 
users  access  to  the  system.  Each  user  generates  and  transmits  the  signals 
independently.  The  transmitter  includes  a  channel  encoder,  a  channel  inter¬ 
leaver  and  an  M-ary  FSK  modulator.  The  channel  code  and  the  interleaver 
are  optional.  For  the  M-ary  FSK  modulation,  all  users  share  the  same  set 
of  frequencies  with  size  M.  Since  there  is  no  cooperation  among  the  T  users 
when  they  generate  and  send  the  signals,  it  is  possible  that  some  of  the  users 
send  the  same  frequency  at  the  same  time,  and  collision  happens.  At  the 
receiver,  the  demodulator  operators  on  the  signals  over  the  M  frequencies. 
If  the  demodulator  outputs  the  number  of  users  that  transmit  on  each  fre- 
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quency  by  detecting  the  energy  level  on  that  frequency  point,  the  channel  is 
called  the  CFSK  channel  with  intensity  information.  In  this  case,  the  output 
of  the  demodulator  is  given  by  a  vector  F,  in  which  the  entries  denote  the 
number  of  transmissions  on  the  frequencies.  On  the  other  hand,  when  the  de¬ 
modulator  only  detects  whether  there  is  any  transmission  on  each  frequency 
point,  the  channel  is  called  the  CFSK  channel  without  intensity  information. 
In  this  case,  the  entry  of  the  vector  F  is  one  if  there  are  transmissions  on 
the  corresponding  frequency,  zero  otherwise.  The  decoder  does  multiuser  de¬ 
tection  on  the  output  of  the  demodulator,  F.  For  each  user,  the  component 
decoder  regards  other  users’  signals  as  useful  information  and  makes  use  of 
them.  The  output  of  the  multiuser  detection  is  the  information  sequence  for 
each  user. 


Figure  2.1:  The  CFSK  system. 
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2.2  The  CFSK  Transmitter 

2.2.1  Overview 

Each  user’s  transmitter  contains  a  channel  encoder,  a  random  channel  inter¬ 
leaver  and  an  M-ary  FSK  modulator.  The  channel  encoder  and  interleaver 
are  optional.  The  code  can  be  block  code  or  convolutional  code.  The  purpose 
for  the  channel  code  is  to  provide  error  control  for  the  system.  Moreover,  the 
application  of  the  channel  code  can  provide  the  user  unique  feature  (UUF) 
[23]  so  that  the  users  can  be  distinguished  at  the  receiver.  When  the  trellis 
code  presents  m  the  system,  the  CFSK  system  can  be  regarded  as  a  trellis 
coded  modulation  system.  The  random  channel  interleaver  after  the  channel 
encoder  is  also  optional.  The  interleaver  is  symbol-based.  It  can  also  provide 
the  UUF.  Purpose  for  the  channel  interleaver  is  to  improve  the  system  perfor¬ 
mance,  which  will  be  discussed  in  Chapter  4.  In  the  M-ary  FSK  modulation, 
the  user’s  coded  sequence  is  mapped  to  the  frequency  signals.  Considering 
that  there  are  M  frequencies,  the  signal  space  is  an  M-dimensional  space, 
where  each  dimension  denotes  a  frequency.  Application  of  different  signal 
mappings  can  also  provide  UUF. 
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2.2.2  User  Unique  Feature  (UUF) 

The  conception  of  UUF  is  given  by  [23].  In  the  multiple  access  system, 
all  users  share  the  common  channel.  Thus,  the  users  need  some  particular 
properties  to  distinguish  themselves  from  others.  For  example,  in  the  SSMA 
systems,  each  user  has  a  specific  spreading  sequence.  In  a  collision  type 
system,  such  as  CFSK,  there  axe  many  ways  to  achieve  UUF.  As  described 
previously,  UUF  can  be  achieved  by  using  different  channel  codes  for  dif¬ 
ferent  users.  Or,  it  can  be  achieved  by  using  different  channel  interleavers 
or  different  signal  mappings  for  different  users.  And,  these  schemes  can  be 
combined  together  to  achieve  better  performance.  Different  UUF  schemes 
require  different  decoder  structures  for  the  multiuser  detection. 

In  this  project,  we  focus  on  the  using  different  channel  codes  for  different 
users  to  achieve  the  UUF.  IVellis  codes,  especially  convolutional  codes,  are 
considered.  All  users  are  using  the  same  signal  mapping  over  a  set  of  M 
frequencies.  Assume  that  (n,  k)  convolutional  codes  are  used  for  all  the  users. 
Then,  we  have  M  =  2”,  i.e.  each  n  bits  of  coded  sequence  are  mapped  to  a 
frequency.  The  channel  interleaver  is  optional  in  this  project,  which  will  be 
applied  in  discussion  of  the  performance  for  the  iterative  detection  in  Chapter 
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2.2.3  Unique  Decodability 


Figure  2.2:  The  super  trellis  code  of  a  T-user  CFSK  system. 


The  transmitter  of  each  user  is  considered  as  a  trellis  coded  modulation 
(TCM)  system,  since  both  block  codes  and  convolutional  codes  are  trellis 
codes.  When  all  users  are  considered  together,  the  transmitter  can  be  re¬ 
garded  as  a  super  trellis  code.  In  Figure  2.2,  the  super  trellis  code  is  de¬ 
scribed.  The  code  is  constructed  by  aligning  all  users’  input  bits  together, 
and  aligning  all  the  memories  together.  When  rate  k/n,  memory  m  convo¬ 
lutional  codes  are  applied,  the  trellis  code  of  the  TCM  has  kT  bits  input 
and  mT  memories  (state  number  is  2-^),  where  T  is  the  number  of  users. 
Since  each  n  bits  of  coded  sequence  are  mapped  to  a  frequency  and  there  are 
M  frequencies,  the  code  symbols  on  the  branch  of  the  super  trellis  can  be 
expressed  by  the  M-ary  vector  F.  The  entry  of  the  F  denotes  the  number 
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of  users  that  transmit  the  frequency.  The  code  symbols  described  by  F 
occupy  the  points  in  an  M-dimensional  signal  space.  For  each  signal  point, 
the  sum  of  its  coordinates  equals  to  the  user  number  T.  As  an  example, 
signal  space  for  a  3-user  3-frequency  system  can  be  described  in  Figure  2.3. 


Figure  2.3:  The  signal  space  for  a  3-user  3-frequency  CFSK  system. 

Unique  decodability  of  the  super  trelUs  code  has  significant  influence  on  the 
performance  of  the  CFSK  system.  The  definition  of  the  uniquely  decodable 
code  (UDC)  for  a  multiple  access  channel  was  first  given  by  Kasami  and  Lin 
in  [24]  in  1976.  Their  discussion  is  based  on  the  block  code  for  the  binary 
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adder  channel.  In  1979,  Chang  and  Weldon  gave  some  construction  methods 
for  the  uniquely  decodable  block  codes  for  the  T-user  adder  channel.  These 
constructions  can  approach  the  maximal  code  rates  given  by  capacity  of  the 
T-user  adder  channel.  Further,  Chang  and  Wolf  put  forward  the  definition 
of  the  uniquely  decodable  block  code  for  the  CFSK  channel  and  gave  the 
construction  methods  for  the  codes,  which  can  approach  the  channel  capacity. 
The  conception  of  the  uniquely  decodable  convolution  code  is  first  given  by 
Peterson  in  1979  in  [25],  where  the  two-user  binary  adder  channel  is  assumed. 
Chevillat  developed  the  results  further  in  [26].  He  gave  the  definition  and 
some  constructions  for  the  uniquely  decodable  trellis  code  for  the  N-user 
adder  channel.  For  the  CFSK  system,  if  every  path  through  the  super-trellis 
has  a  distinct  code  sequence,  then  the  code  set  is  uniquely  decodable.  In 
this  situation,  each  code  path  has  to  start  from  the  state  zero  and  end  at  the 
state  zero.  Each  codeword  in  the  code  sequence  is  given  by  F.  Brute  force 
checking  shows  that  the  two  codes  given  in  Figure  2.4  is  uniquely  decodable. 


01 


Figure  2.4:  A  two-user  super-trellis. 
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It  has  been  shown  that  to  construct  a  set  of  uniquely  decodable  codes  with 
error  control  for  a  large  number  of  users  is  extremely  difficult  [27]  [28].  For 
later  discussion,  the  recursive  systematic  convolutional  codes  (RSC)  are  gen¬ 
erated  randomly  for  the  users.  An  important  reason  for  applying  RSC  is 
that  the  systematic  convolutional  codes  with  canonical  controller  form  are 
proved  to  be  non-catastrophic  [29]  [30]  .  Since  RSC  is  of  such  form,  the  ran¬ 
domly  generated  codes  axe  non-catastrophic  automatically.  Also,  a  couple  of 
generation  rules  [22]  can  guarantee  the  uniform  implementation  for  all  users. 

In  Figure  2.5,  the  code  set  for  the  3-user  system  is  not  uniquely  decodable. 
For  instance,  the  three  state  sequences  through  the  trellises  are 

code{a)  :  00, 00, 10, 01, 00, 00, 10, 11, 01, 00, 
code{b) :  00, 00, 00, 10, 11, 01, 10, 11, 01, 00, 
code{c) :  00, 10, 01, 10, 11, 01, 10, 01, 00, 00.  (2.1) 

The  corresponding  output  sequence  on  the  super-trellis  is 

(2, 0, 1, 0),  (1, 0, 1, 1),  (1, 1, 0, 1),  (1, 1, 1, 0),  (1, 2, 0, 0), 

(1,1, 0,1),  (0,1, 1,1),  (0,1, 1,1),  (1,0, 1,1). 


(2.2) 
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There  is  another  set  of  state  sequences  through  the  trellises: 

codeia)  :  00, 00, 10, 11, 11, 11, 11,  H,  01, 00, 
code{b) :  00, 00, 00, 00, 00, 00, 10, 01, 00, 00, 
codeic)  :  00, 10, 01, 10, 01, 10, 11, 11, 01, 00,  (2.3) 

with  the  same  output  sequence  given  in  (2.2).  Hence,  whenever  the  code 
sequence  in  (2.2)  is  received,  one  of  paths  in  (2.1)  and  (2.3)  has  to  be  decided 
randomly.  So,  it  is  always  possible  that  decoding  error  will  happen.  And, 

this  will  lead  to  an  error  floor  for  the  bit  error  performance  of  the  CFSK 
system. 


Figure  2.5:  An  code  set  for  the  3-user  CFSK  system,  which  is  not  uniquely 
decodable. 


The  effects  of  the  non-uniquely  decodable  codes  will  be  discussed  in  Chapter 
5.  In  Chapter  6,  some  constructions  of  uniquely  decodable  codes  will  be 


discussed. 
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2.3  The  CFSK  Channel 

2.3.1  Overview 

For  the  M-aiy  FSK  modulation,  let  the  set  of  M  frequencies  be 

=  {/o,  /i,  •  •  • ,  fu-i}-  (2.4) 

Then,  the  output  signal  from  the  channel  has  the  form 

M-l 

sm  = 

where  Hj  is  the  multiplicity  factor  that  denotes  the  number  of  users  that 
transmit  frequency  E  is  the  signal  energy.  is  a  function  of  fXj  that 
depends  on  the  synchronization  of  the  users,  which  will  be  discussed  later. 
6j  is  a  uniformly  distributed  the  phase  for  frequency  fj.  nj  is  the  additive 
noise.  In  Chapter  3,  specific  math  model  for  the  CFSK  will  be  given  for  the 
synchronous  and  asynchronous  channels. 

2.3.2  Synchronization 

In  Figure  2.6,  three  synchronization  modes  are  described.  For  the  full  syn¬ 
chronous  mode,  all  users  start  transmission  at  the  same  time  and  with  the 
same  phase  for  each  frequency.  For  the  “asynchronous  over  time  step”  mode, 
the  users  may  start  transmission  at  different  time.  But,  they  transmit  a 


(l>{fjtj)y/Ecos{2nfjt  +  9j)  +  nj  ,  (2.5) 
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Figure  2.6:  Synchronization  modes  on  the  CFSK  channel. 
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frequency  with  the  same  phase.  Thus,  the  symbols  are  synchronized.  If  ad¬ 
ditional  communication  protocol  can  guarantee  that  the  receiver  knows  when 
each  user  starts  transmission,  then  the  above  two  modes  axe  the  same  for  the 
receiver.  Thus,  an  assumption  for  later  discussion  is  that  the  the  above  two 
modes  are  all  synchronous  case.  The  third  case  in  Figure  2.6  is  the  asyn¬ 
chronous  mode.  In  this  mode,  the  users  transmit  their  symbols  with  different 
phases  no  matter  they  start  transmission  together  or  not.  In  (2.5),  the  func¬ 
tion  ^{Hj)  depends  on  the  synchronization  mode.  For  the  synchronous  mode, 
=  tij-  This  is  because  that  all  users  transmit  fj  with  the  same  phase, 
which  is  equal  to  6j,  thus  the  result  signal  on  fj  is  just  a  sum  of  all  the  sig¬ 
nals  of  the  users.  For  the  asynchronous  case,  the  value  of  (j>{iij)  is  a  random 
variable  with  variance  HjE/2,  which  will  be  discussed  in  Chapter  3. 

2.4  CFSK  Receiver  Structures 

2.4.1  Overview 

The  frequency  spectrum  of  the  channel  output  signals  only  occupies  on  single 
frequency  points.  Thus,  a  simple  demodulation  method  can  be  achieved  by 
using  the  discrete  Fourier  transformation  (DFT).  Since  M  frequencies  are 
used  in  the  FSK  modulation,  the  signal  has  a  bandwidth  of  M.  So,  2M 
points  need  to  be  sampled  in  each  symbol  period.  The  sampled  points  are 
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feed  into  the  2M-point  DFT.  Then,  the  energy  level  and  phase  for  each 
frequency  are  output  from  the  DFT  function.  For  frequency  suppose  that 
(2.5)  gives  the  channel  signal.  Then,  the  DFT  output  on  fj  can  be  given  by 
a  phasor  expression 

=  Pj^^j  =  [<f>{fij)fEiej.  (2.6) 

For  synchronous  case,  pj  equals  to  p)E  without  considering  the  noise.  For 
asynchronous  case,  pj  is  a  random  variable  with  expectation  pjE.  Thus,  the 
output  of  DFT  over  all  frequencies  can  be  expressed  by  the  vector 

^  =  {Po,Pi,.. .  ,Pm-i)-  (2.7) 

The  multiuser  detection  takes  the  vector  F  as  input.  For  each  user,  the 
information  from  other  users  is  not  treated  as  interference.  The  information 
from  other  users  are  applied  in  decoding  each  user’s  information.  Generally, 
there  are  two  specific  approaches  for  multiuser  detection. 

2.4.2  Optimal  Detection  of  CFSK 

The  optimal  detection  gives  the  best  performance  for  the  multiuser  detection 
for  the  CFSK  system.  It  operates  on  the  super  trellis  of  the  system.  By 
maximum  likelihood  sequence  detection  (MLSD)  on  the  super  trellis  [14], 
the  decoder  can  provide  the  best  estimation  for  the  code  path  through  the 
super  trelhs.  From  a  specific  state  path  through  the  super  trellis,  each  user’s 
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state  path  through  his  trelHs  can  be  uniquely  decided.  Thus,  each  user’s 
coded  sequence  and  information  sequence  can  be  decided.  Note  that  when 
the  code  set  is  not  uniquely  decodable,  there  will  be  two  or  more  paths  on 
the  super  trellis  with  the  same  code  sequence.  In  such  circumstance,  the 
decoding  algorithm  has  to  choose  a  path  randomly.  So,  it  is  always  possible 
that  there  are  incorrect  estimations  of  the  state  path  without  considering  the 
noise  on  the  channel. 

Though  MLSD  approach  can  provide  the  best  detection  results,  its  compu¬ 
tational  complexity  is  not  polynomial  (NP-complexity),  which  means  that 
the  computational  complexity  of  MLSD  will  increase  with  the  user  number 
exponentially.  The  reason  is  that  the  super  trellis  aligns  the  users’  input 
and  memory  together.  Thus,  the  possible  size  of  input  alphabet  and  state 
number  will  increase  exponentially  when  the  user  number  increases.  As  a  re¬ 
sult,  iterative  detections  are  of  interests,  which  have  polynomial  computation 
complexity  and  can  approach  the  performance  of  MLSD. 

2.4.3  Iterative  Detection  of  CFSK 

Unlike  the  optimal  detection,  which  operates  on  the  super  trellis  and  has  a 
single  decoder,  the  iterative  detection  approach  decodes  each  user’s  informa¬ 
tion  by  their  component  decoders.  After  each  iteration,  the  users  exchange 
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their  updated  information  with  each  others.  The  obtained  information  from 
others  are  applied  in  the  following  iteration.  After  enough  iterations,  the  de¬ 
coding  results  for  all  the  users  are  supposed  to  approach  the  optimal  estima¬ 
tion.  Two  known  iterative  detection  approaches  are:  interference  cancellation 
and  consensus  decoder,  which  treat  the  updated  information  differently. 


2.4.4  Interference  Cancellation  (IC) 


Figure  2.7:  The  interference  cancellation  decoder. 


The  structure  of  interference  cancellation  is  given  in  Figure  2.7.  In  each  it¬ 
eration,  the  component  decoders  update  the  estimation  of  the  transmitted 
frequencies.  At  the  input  of  the  component  decoders,  the  estimations  from 
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other  users  axe  subtracted  from  the  vector  F,  so  that  the  input  of  the  com¬ 
ponent  decoder  is  an  updated  channel  signal  for  a  specific  user.  Information 
form  other  users  is  canceled  from  the  channel  signal  F  as  interference.  And, 
each  component  decoder  treats  the  channel  as  a  single  user  channel  by  using 
single  user  channel  metric.  This  make  it  seems  that  IC  is  a  kind  of  single 
user  detection.  However,  for  eaeh  user,  the  single  user  channel  is  obtained 
by  using  the  updated  estimations  of  the  frequencies  from  other  users.  In  this 
sense,  IC  is  regarded  as  an  iterative  multiuser  detection. 

In  [23],  performance  of  IC  for  TCMA  system  is  given,  which  shows  that 
the  performance  for  IC  is  much  worse  than  the  optimal  performance.  This  is 
because  that  the  component  decoders  treat  the  channel  as  single  user  channel 
and  use  single  user  channel  metric.  So,  this  project  will  not  put  much  energy 
on  the  IC  approach. 

2.4.5  The  Consensus  Decoder 

The  structure  of  consensus  decoder  is  given  by  [18].  In  Figinre  2.8,  the  ba¬ 
sic  structure  of  consensus  decoder  is  described.  The  consensus  decoder  is 
composed  of  the  users’  component  decoders.  The  component  decoder  takes 
the  DFT  output  F  as  input  and  decode  each  user’s  information  separately. 
Each  component  decoder  outputs  the  estimations  for  each  user’s  symbols  and 
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Figure  2.8:  The  consensus  decoder. 
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the  a  posteriori  information  on  the  symbols.  For  instance,  is  the  a 

posteriori  probabihty  for  frequency  output  by  user  The  a  posteriori 
information  is  feed  to  other  users’  component  decoder  as  a  priori  information. 
The  a  priori  information  from  other  users  are  applied  to  update  the  channel 
metric  during  each  user’s  decoding  procedure.  After  enough  iterations,  the 
estimation  results  for  each  user  is  supposed  to  approach  the  optimal  results. 
In  Chapter  4  the  consensus  decoders  using  MAP  decoder  and  Turbo  itera¬ 
tive  decoder  are  introduced.  The  performance  for  the  consensus  decoders  is 
discussed  in  Chapter  5. 


Chapter  3 

Theory  and  Performance  of  the 
CFSK  Channel 

3.1  Introduction 

The  collision  frequency  shift-keyed  (CFSK)  channel  is  based  on  the  channel 
model  described  by  Chang  and  Wolf  as  the  T-user  M-frequency  multiple- 
access  channel  with  and  without  intensity  information  [17].  They  gave  out 
the  capacities  of  the  noiseless  channels  with  and  without  intensity  informa¬ 
tion.  A  formal  description  of  the  CFSK  channel  was  given  by  Grant  and 
Schlegel  in  [21]  and  [18].  Starting  from  results  in  [17],  Grant  derived  the 
approximated  capacity  for  the  noiseless  CFSK  channel  without  intensity  in¬ 
formation  in  [21],  which  approaches  Chang  and  Wolf’s  result  asymptotically. 
Furthermore,  Grant  provided  the  upper  and  lower  bounds  of  the  CFSK  chan¬ 
nel  without  intensity  information  in  the  presence  of  noise  and  undetected  in- 
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terference.  These  results  set  up  a  profound  theoretical  basis  for  the  design  of 
a  multiple-access  communication  system  based  on  the  CFSK  channel.  Due 
to  the  complexity  of  the  CFSK  channel  with  intensity  information,  Grant 
did  not  give  out  the  capacity  of  the  noisy  CFSK  channel  with  intensity  in¬ 
formation.  Moreover,  all  the  former  capacity  results  are  assuming  that  the 
users  that  access  to  the  channel  send  their  symbols  synchronously,  i.e.  the 
symbols  on  the  channel  are  aligned  with  each  other  perfectly.  In  practical  ap¬ 
plications,  the  synchronization  is  hard  to  be  achieved  among  multiple  users. 
For  an  asynchronous  CFSK  system,  users  send  their  symbols  without  any 
cooperation,  which  means  that  the  phase  of  the  FSK  signal  sent  by  a  user 
on  a  particular  frequency  is  different  from  those  sent  by  other  users.  This 
is  similar  to  introducing  a  noise  to  the  synchronous  channel.  Thus,  different 
capacity  results  of  the  CFSK  channel  are  anticipated. 

In  [22]  and  [31],  mathematical  models  for  the  synchronous  and  asynchronous 
CFSK  channels  are  set  up.  These  models  provide  a  new  approach  to  research 
the  asynchronous  CFSK  channel.  On  the  other  hand,  the  development  of 
computer  technology  make  it  possible  to  study  the  more  complicated  chan¬ 
nels,  such  as  the  CFSK  channel  with  intensity  information.  In  this  work, 
efforts  are  made  to  study  the  sum  capacities  for  the  CFSK  channels  with  and 
without  intensity  information  in  noisy  or  asynchronous  scenarios.  Inspired 
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by  ideas  in  [21],  analytical  models  are  set  up  so  that  the  bounds  of  the  sum 
capacities  are  obtained.  The  upper  and  lower  bounds  of  the  sum  capacity 
of  the  noisy  synchronous  CFSK  channel  with  intensity  information  are  com¬ 
puted  in  section  3,  which  provides  a  computer  approach  for  computing  the 
sum  capacity.  Then,  the  upper  bounds  of  the  noiseless  asynchronous  CFSK 
channel  with  intensity  information  is  obtained  by  the  computer  approach. 

By  the  analytical  methods,  an  upper  bound  for  the  noiseless  asynchronous 
CFSK  channel  without  intensity  information  is  derived.  Results  show  that  all 
the  bounds  are  tight  enough  for  large  user  numbers.  These  capacity  results 
are  compared  together,  which  provides  significant  directions  in  designing  a 
CFSK  system  for  practical  applications. 

3.2  The  CFSK  Channel  Distribution 

3.2.1  Distribution  over  the  Synchronous  Channel 

According  to  the  channel  model  in  [22]  and  [31],  the  output  signal  of  the 
noisy  synchronous  CFSK  channel  on  frequency  j  can  be  described  by 

Xj  =  fXjy/Ecos{(f)j)  -f  n^, 

Vj  =  fijVEsin{<f)j)  +  Uy,  (3.1) 

Pj  =  x]  +  y], 

where  fXj  is  the  number  of  users  that  transmit  frequency  j,  E  is  the  symbol 
energy,  (f>j  is  the  phase  with  which  all  the  users  transmit  the  signal,  and  Uy 
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are  the  the  Gaussian  noise  with  double-sided  power  spectral  density  (P.S.D.) 

Thus,  pj  is  the  detected  energy  on  the  output  of  the  channel  at  frequency 
j.  The  distribution  of  pj  is  given  by  a  non-central  chi-square  distribution, 
which  is  further  approximated  by 


{pjVE-^y 

N„ 


(3.2) 


The  relationship  between  the  above  fpj{pj)  and  is  described  in  Figure  3.1 
for  =  {0, 1, . . . ,  8}  when  the  signal-to-noise  ratio  (SNR)  is  8dB. 


Probability  Density  Function  for  The  Symbol  Energy  of  The 
Noisy  Synchronous  CFSK  Channel  at  8dB 


Figure  3.1:  P.S.D.  of  symbol  energy  of  the  noisy  synchronous  CFSK  channel. 


CHAPTER  3.  THEORY  AND  PERFORMANCE  OF  THE  CFSK  CHANNEL32 

3.2.2  Distribution  over  the  Asynchronous  Channel 


For  the  asynchronous  CFSK  channel  without  considering  the  noise,  the  out¬ 
put  signal  given  by  [22]  and  [31]  is 

xj  =  VEY:^l,cos{cl>j,i), 

yj  =  (3.3) 

pj  = 

where  denotes  the  phase  with  which  the  user  i  transmits  frequency  j. 
Since  all  users  transmit  the  signals  without  any  cooperation,  the  distribution 
of  the  phase  can  be  assumed  uniform  between  0  and  27r  as 


^  0  <  4>j,i  <  27r 

0  otherwise 


(3.4) 


Theorem  1  (Variance  of  noiseless  asynchronous  signals):  The  variance  of  Xj 

and  yj  given  m  (3.3)  is  and  the  signal  can  be  approximated  by  a  noisy 
synchronous  signal  as 

Xj  =  Pjy^cos{<f>j)+neff,^, 

Vi  =  yj^sin{(f)j)  +  neff,y,  (3.5) 

Pi  =  a;J  +  y?, 

where  the  and  are  the  effective  Gaussian  noise  with  double-sided 

P.S.D. 

The  proof  of  Theorem  1  is  given  in  Appendix  I.  From  Theorem  1,  the  distri¬ 
bution  oipj  over  the  noiseless  asynchronous  CFSK  channel  can  be  obtained 
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in  a  similax  way  as  the  one  over  the  noisy  synchronous  CFSK  channel.  Thus, 


UAp,)  =  - 


HjE 


(3.6) 


The  relationship  between  fpj  (pj)  and  is  described  in  Figure  3.2  for  pj  = 


Probability  Densi^  Function  for  The  Symbol  Energy  of  The 
Noiseless  Asynchronous  CFSK  Channel 


=  (x^  +  y®)’®  (p;  tfie  symbol  energy) 


Figure  3.2:  P.S.D.  of  symbol  energy  of  the  noiseless  asynchronous  CFSK 
channel. 


In  both  Figure  3.1  and  Figure  3.2,  it  can  be  observed  that,  the  curves  of  the 
distributions  extend  to  the  area  where  y/p  is  less  than  0.  This  is  the  result  of 
the  approximation  of  the  non-central  chi-square  distribution,  which  will  not 
affect  later  discussion  since  all  integrations  are  from  —  (X)  to  positive  values  of 
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y/p.  The  distribution  when  fXj  is  1  is  not  given  in  Figure  3.2,  for  the  reason 
that  there  is  no  collision  on  frequency  j  when  pj  =  1,  and  distribution  of 
the  asynchronous  channel  is  the  same  as  that  of  the  synchronous  channel.  It 
can  also  be  observed  that,  for  large  pj  in  Figure  3.2,  the  distributions  overlap 
with  each  other  seriously.  This  implies  that  it  is  much  less  precise  to  estimate 
energy  on  the  asynchronous  channel.  However,  it  would  be  better  to  detect 
whether  there  is  transmission  on  the  frequency  or  not.  So,  it  is  reasonable 
to  anticipate  that  the  performance  of  the  asynchronous  CFSK  channel  with 
intensity  information  may  be  worse  than  that  of  the  asynchronous  CFSK 
channel  without  intensity  information  when  the  system  is  over-loaded.  Later 
results  will  prove  this  anticipation. 


3.3  Performance  of  the  Synchronous  CFSK 
Channel 

[17]  provided  the  sum  capacities  of  the  noiseless  synchronous  CFSK  channel 
with  and  without  intensity  information.  These  results  can  be  observed  in 
Figure  1.1  when  they  are  compared  Avith  the  sum  capacities  of  the  CDMA 
systems  [18].  Grant  provided  a  tight  approximation  of  the  capacity  for  the 
noiseless  synchronous  CFSK  channel  without  intensity  information  in  [21]. 
Furthermore,  upper  and  lower  bounds  for  the  sum  capacity  of  the  noisy 
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synchronous  CFSK  channel  without  intensity  information  is  also  given  in 
[21].  In  this  context,  discussion  on  the  bounds  for  the  noisy  synchronous 
CFSK  channel  with  intensity  information  is  provided. 

Theoretically,  the  size  of  the  output  alphabet  from  the  noisy  CFSK  channel 
with  intensity  information  is  infinite.  This  makes  it  extremely  difficult  to 
compute  the  exact  capacity  of  the  noisy  CFSK  channel  with  intensity  infor¬ 
mation.  Grant’s  analytical  channel  model  in  [21]  provides  an  effective  way 
for  computing  the  upper  and  lower  bounds  of  the  noisy  CFSK  channel.  The 
idea  of  the  model  is  to  decompose  the  noisy  channel  into  several  concatenated 
sub-channels.  Each  sub-channel  generates  specific  errors  that  would  appear 
on  the  original  channel.  The  model  is  modified  in  this  context  for  computing 
the  upper  and  lower  bounds  of  the  noisy  synchronous  channel  with  intensity 
information. 

An  analytical  model  for  the  noisy  synchronous  channel  with  intensity  infor¬ 
mation  is  given  in  Figure  3.3.  X,  Y,  U,  V  and  W  are  the  input  and  output 
alphabets  of  the  sub-channels.  Each  element  of  the  alphabets  is  a  1  x  M  vec¬ 
tor  in  which  each  entry  denotes  the  number  of  transmissions  (occupations) 
on  the  corresponding  frequency.  For  all  x  e  X,  let  x  be 


X  (Xi ,  X2 ,  •  •  •  ,  ) . 


(3.7) 
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And,  for  all  ueU,  let  u  be 


^  —  (^11^2,  .  .  .  iUm). 


(3.8) 


Then,  the  single-hop  denotes  the  kind  of  error  that,  when  the  Xj  users  trans¬ 
mit  on  frequency  j,  the  output  of  the  single-hop  channel  on  frequency  j 
satisfies  that 


/  ^3  no  error  happens 
^  ^  Xj  ±  1  otherwise  (3-^) 

The  alphabet  V,  which  is  the  output  of  the  noisy  channel,  can  include  other 

errors  defined  by  n-hop  for  n  >  1.  Let  u  €  ^  be  defined  by 


(3.10) 


Then,  the  n-hop  error  happens  when 


Vj  =  Xj  db  n. 


(3.11) 


In  Figure  3.3,  the  “symmetrifying  channel”  takes  the  noised  signal  as  input. 
At  the  output  of  this  channel,  all  kinds  of  errors  will  have  the  same  probability 


Figure  3.3:  The  analytical  channel  model  for  the  CFSK  channel  with  inten¬ 
sity  information. 
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The  probability  of  single-hop,  denoted  by  q  ,  can  be  computed  by  the  distri¬ 
bution  given  by  (3.2).  Then,  let  the  probability  Psym  be  equal  to  q.  According 
to  the  data  processing  inequality  [32],  the  sum  capacity  of  the  noisy  channel 
has  the  following  upper  and  lower  bounds: 

I{W\X)  <  Cnaisy,sum  =  IiV\X)  <  I{U\X).  (3.12) 

A  lower  bound  can  be  computed  from  I{W\X)  as  given  in  the  following 
Theorem.  A  proof  of  the  similar  problem  can  be  found  in  [21]. 

Theorem  2:  a  lower  bound  of  the  sum  capacity  of  the  noisy  synchronous 
CFSK  channel  with  intensity  information  is  given  by 

Cnoisy,sum  =  I{V\X)  >  I{W\X)  >  H{Y)  -  Mh{q)  >  Csum  -  Mh{q),  (3.13) 

where  Cs^m  is  the  sum  capacity  of  the  noiseless  synchronous  CFSK  channel 
with  intensity  information,  M  denotes  the  number  of  frequencies  used  in  the 
system,  h{q)  is  the  binary  entropy  function  given  by 

Hq)  =  -qiog{q)  -  (1  -  q)log{l  -  q).  (3.14) 

The  value  I{U\X)  provides  an  upper  bound  for  the  sum  channel  capacity, 
which  can  be  computed  by 


I{U\X)  =  H{U)  -  H{U\X). 


(3,15) 
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The  entropy  H{U)  is  given  by 


^  -pu{u)log{pu{u)).  (3.16) 

ueu 

The  mutual  information  is  given  by 


H{U\X)  —  ^  Pu,x{u,x)I{u\x) 

ueu,xex 


Considering  that 


^  -Px{x)pu\x{u\x)log{pu\x{u\x)). 

ueu.xex 

(3.17) 


Pu{u)  -  Y^Px{x)pu\xiu\x),  (3.18) 

xex 

it  is  important  to  compute  the  channel  transition  probability  Pu\x{u\x).  Since 
U  IS  the  output  of  the  sub-channel  with  single-hop  error,  pu\x{u\x)  is  equal 
to  0  when  an  output  u  cannot  be  generated  by  a  given  input  x.  Thus,  for  a 
certain  output  u  of  the  sub-channel,  a  set  can  be  constructed  as 


S  {x\x  €  X,  u  is  the  output  of  the  sub-channel  by  input  x.}.  (3.19) 

As  a  result,  (3.18)  becomes 


Pu{u)  -  y^^Px{x)pu\x(u\x).  (3.20) 

xes 

For  all  X  G  5,  the  probability  pu\x{u\x)  is  given  by 


PU\x{tJ'\x)  =  g'^(“>®)(l  _ 


(3.21) 


where  the  function  W{u,x)  denotes  the  number  of  positions  in  u  that  are 
different  from  those  in  x.  Assume  that  all  users  send  the  frequencies  with 
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equal  probability,  which  is  Then,  the  probability  of  input  sjonbol  x  is 
computed  as 


px{x) 


MM 


(3.22) 


By  computer  approach,  it  is  not  difficult  to  construct  the  set  of  all  input  sym¬ 
bols  X.  By  the  rule  of  single-hop,  the  set  U  can  also  be  constructed.  When 
the  frequency  number,  M,  and  the  user  number,  T,  are  not  very  large,  the 
time  spent  on  these  constructions  is  not  huge.  Following  the  computations 
given  by  (3.15)-(3.22),  the  upper  bound  can  be  obtained.  In  Figure  3.4,  the 
upper  bounds  and  lower  bounds  for  the  system  when  SNR  is  8dB  are  given. 
It  can  be  found  that  the  bounds  are  tighter  when  the  number  of  users  gets 
large.  When  SNR  increase,  the  bounds  will  also  become  tighter.  This  can 
be  observed  in  Figure  3.5,  where  the  bounds  are  given  when  SNR  is  12dB. 
The  sum  capacities  of  the  noiseless  channels  [17]  is  also  given  in  Figure  3.4 
and  3.5.  For  high  SNR’s,  those  bounds  approach  the  sum  capacity  of  the 


noiseless  channel. 
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bounds  for  capacity  of  noisy  synchronous  CFSK 
channel  with  intensity  information  at  8dB. 
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Upper  arxJ  Lower  Bounds  for  Capacity  of  Noisy  Synchronous 
CFSK  Chanel  with  Intensity  Infomnation  at  12dB 


Figure  3.5:  Upper  and  lower  bounds  for  capacity  of  noisy  synchronous  CFSK 
channel  with  intensity  information  at  12dB. 
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3.4  Performance  of  the  Asynchronous  CFSK 
Channel 

3.4.1  The  Noiseless  Asynchronous  Channel  with  In¬ 
tensity  Information 


To  analyze  the  capacity  of  the  noiseless  asynchronous  CFSK  channel,  the 
analytical  model  given  in  Figure  3.3  is  still  applied.  The  basic  idea  of  com¬ 
puting  the  upper  bound  of  the  capacity  follows  the  one  given  in  the  previous 
section.  The  only  difference  exists  in  the  probability  pu\x{u\x),  since  the 
distribution  of  the  channel  is  given  by  (3.6).  From  (3.6),  a  transmission 
probability  matrix  can  be  constructed  as 


P(0|0)  p(l,0)  ...  p{T  +  j3\{))  \ 

p(0|l)  p(l,l)  ...  p{T  +  p\l) 

Vp(0|T)  p(l,r)  ...  p{T  +  l3\T)  } 


(3.23) 


where  p(wj  |a:j)  denotes  that  probability  that  Uj  transmissions  are  detected  at 


the  output  of  the  single-hop  channel  given  that  Xj  users  transmit  on  frequency 
j.  /?  is  a  constant,  which  is  large  enough  that  p{T  +  0\xj)  can  be  ignored  for 
all  0  <  oTj  <  T.  Then,  the  probability  Pu\x{'u\x)  is  given  by 


M 

Pu\x{u\x)  =  JJp(Uj|Xj). 

i=i 


(3.24) 


By  computer  approach,  the  alphabets  X  and  U  are  constructed,  and  com¬ 
putations  given  by  (3.15)-(3.20),  (3.24)  and  (3.22)  are  followed.  The  upper 
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Upper  Bound  for  Capacity  of  Noiseless  Asynchronous 
CFSK  Channel  with  Intensity  Information 


Figure  3.6:  Upper  bounds  for  capacity  of  noiseless  asynchronous  CFSK  chan¬ 
nel  with  intensity  information. 


CHAPTER  3.  THEORY  AND  PERFORMANCE  OF  THE  CFSK  CHANNEIA4 

bounds  when  2  <  M  <  8  are  described  in  Figure  3.6.  It  can  be  observed 
that  when  user  number,  T,  becomes  large,  especially  when  T  >  M,  the  up¬ 
per  bound  of  the  capacity  decreases  drastically.  This  is  compared  to  the 
capacity  of  the  noisy  synchronous  channel,  which  keeps  increasing  when  user 
number  increases.  And,  the  upper  bound  tends  to  0  when  user  number  is  ex¬ 
tremely  large.  So,  the  asynchronous  channel  is  unusable  under  heavy-loaded 
conditions. 

Since  the  upper  bound  described  in  Figure  3.6  contains  enough  information 
to  show  the  performance  of  the  channel,  and  it  is  too  complicated  to  compute 
the  lower  bound  from  the  distribution  given  by  (3.6),  the  lower  bound  of  the 

noiseless  asynchronous  channel  with  intensity  information  is  not  given  in  this 
context. 

3.4.2  The  Noiseless  Asynchronous  Channel  without  In¬ 
tensity  Information 

To  compute  the  capacity  of  the  noiseless  asynchronous  CFSK  channel  with¬ 
out  intensity  information,  an  analytical  model  is  given  in  Figure  3.7.  By 
this  model,  only  upper  bound  of  the  capacity  is  available,  which  is  given 
by  I{U\X).  On  the  first  sub-channel,  errors  only  happen  when  there  are  2 
transmissions  on  frequency  j  {xj  =  2).  Since  the  channel  is  without  intensity 
information,  the  errors  are  always  2-hop,  i.e.  a  jump  from  double-occupation 
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The  noiseless  asynchronous  CFSK  channel  without  intensity  information 


Figure  3.7:  The  analytical  channel  model  for  the  noiseless  asynchronous 
CFSK  channel  without  intensity  information. 


Before  computing  the  upper  bound  of  the  channel,  some  variables  are  defined 
as 

m:  random  variable  counts  active  frequencies  at  the  output  U ; 

V.  random  variable  counts  active  frequencies  at  the  input  X\ 
s:  random  variable  counts  single  occupied  frequencies  at  the  input  X\ 
d:  random  variable  counts  double  occupied  frequencies  at  the  input  X. 

An  upper  bound  is  given  by  the  following  theorem,  which  is  proved  in  Ap¬ 
pendix  II. 


Theorem  3:  The  upper  bound  of  the  noiseless  asynchronous  channel  without 
intensity  information  is  given  by 


^  \  ^  T 

I{U\X)  =  Yj>M{m)log (  I  +  HQ)'^d'^Piii)PD\iid\i),  (3.25) 

where 

T 

T  f 

PM{m)  =  Y^^Piii)pD\i{s\i)pM\D,iim\d,  i), 

i—m  rf=l 


(3.26) 
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PD\i{d\i)  = 


EU  C)  Os'  t/)  nfrj 


(3.27) 

(3.28) 


q  IS  the  probability  of  double-occupation  to  inactive  error,  which  is  equal  to 
p(0|2)  in  (3.23).  4^  denotes  that  number  of  ways  that  put  k  balls  into  n 
bins  such  that  every  bin  contains  at  least  t  balls,  is  computed  in  an 
iterative  way  as 


(3.29) 

(3.30) 

(3.31) 

(3.32) 


The  upper  bounds  are  described  in  Figure  3.8  for  2  >  M  >  16.  By  com¬ 
paring  with  the  capacity  of  the  noiseless  synchronous  CFSK  channel  without 
intensity  information,  it  is  found  that  the  capacity  of  the  asynchronous  chan¬ 
nel  decreases  when  the  number  of  users  gets  large,  especially  when  T  >  M. 
For  small  T,  the  capacity  of  the  asynchronous  channel  decreases  a  little.  By 
results  m  Figure  3.6  and  3.7,  it  is  found  that  the  upper  bounds  of  the  chan¬ 
nel  without  intensity  information  is  generally  higher  than  the  bounds  of  the 
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channel  with  intensity  information.  This  is  reasonable,  since  it  is  less  precise 
to  estimate  the  energy  level  of  the  frequency  than  to  detect  the  active  or 
inactive  of  the  frequency  when  the  users  are  not  synchronized. 


Figure  3.8:  Upper  bounds  for  capacity  of  noiseless  asynchronous  CFSK  chan¬ 
nel  without  intensity  information. 


3.5  Summary 

From  the  upper  bounds  given  in  previous  sections,  it  can  be  concluded  that 
the  asynchronous  CFSK  channels  have  poor  performance  under  heavy-loaded 
conditions.  For  the  channel  with  intensity  information,  the  capacity  ap¬ 
proaches  0  when  user  number  is  much  larger  than  the  frequency  number.  It 
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is  even  worse  than  the  channel  without  intensity  information,  since  the  dis¬ 
tributions  of  signals  with  large  occupation  numbers  overlap  with  each  other 
seriously.  In  this  case,  it  is  better  to  detect  the  active  or  inactive  of  each 
frequency.  When  user  number  is  smaller  than  the  frequency  number,  the 
synchronous  channels  and  the  asynchronous  channels  behave  similarly.  Ap¬ 
plication  of  the  asynchronous  channels  when  the  spreading  factor  (given  by 

f )  IS  between  1  and  10  is  much  more  significant,  since  the  channels  are  more 
efficient  under  such  conditions. 


3.6  Appendix  I:  Proof  of  Theorem  1 


Consider  the  variable  y  which  is  a  function  of  the  phase 

y  =  Vlsm(0),  (3  33) 

where  the  probability  density  function  (P.D.F.)  of  the  phase  $  is  given  by 
Ui(f>)  in  (3.4).  When  y  <  -y/E,  the  probability  function  of  y  is 


■^y(y)  =  0. 


When  y/E  <  y  <  0,  the  probability  function  of  y  is 


Fviy) 


J‘27r  “  arcsin  (y  /  y/E) 
TT— arcsm(y/v/S) 
7r+2arcsm(y/\/^^ 
27r 


(3.34) 


(3.35) 
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When  0  <y  <  y/E, 


Fyiy)  -  fo  /-rarcsm(j//v/E)+7r/*(*^)‘^^+ 


_  arcsin(i//y/E)  ,  arcsin(p/\/B)  ,  1 

27r  ^  2n  “*2 

_  7r-l-2arcsin(y  /  vE) 

27r 


When  y  >  y/E, 


Pyiv)  =  1- 


To  sum  up  the  above  results,  the  probability  function  of  y  is 

fviy)  =  y=\ 


0  y  <  —y/E 

n+2arcsin{y/ym  ^  ^ 

1  y>  y/E 


Thus,  the  P.D.F.  of  Y  from  Fyiy)  is  computed  as 

fy(y)  =  ^ 


~y/E  <y  <  y/E 


0  otherwise 

Now,  introduce  the  moment  generation  function  of  T  as 


^y(s)  =  E[e^‘^] 


=  I^e"^fYiy)dy 

=  II 


e^y 


;dy- 


(3.36) 


(3.37) 


(3.38) 


(3.39) 


(3.40) 


nVE^/l-y^/E 

Let  the  variable  'if  be  the  sum  of  n  variables  with  P.D.F.  fviy)-  Then,  we 
have 

n  n 

if  =  '^yi  =  '^  y/Esin{<f>i),  (3.41) 


1=1  i=l 
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where  the  P.D.F.  of  is  /$((/>)  in  (3.4).  So,  the  moment  generation  function 
of  ^  can  be  obtained  from  $y(s)  as 

rVE 


^^(s)  =  [$y(s)r  = 


/y/t; 

Ve  tta 


^sy 


\dy 


(3.42) 


-Ve  T^y/Eyjl  -  y’^/E 
From  the  moment  generation  function  $4,(5),  the  expectation  and  variance 
of  ^  can  by  obtained.  Firstly,  the  expectation  of  ^  is 


m  =  ^ 


5=0 


=  n 


_ _ 

T^VE^l-y^lE 


1  n~l 


dy 


VVE%-y^lE^'^ 


_ 1 

nVEy/r^^/E 


T  n— 1 


=  n 

=  0. 

Then,  the  variance  of  ^  is 


dy 


s=0 


(3.43) 


uar[^]  =  £^[(^  -  £;[^])2] 


as2 


=  n(n-l)  r^..  dv  r%/g  ye>v 

^VEy/i-yyE^y\  [J-Ve~^^i_^2/e 

r  77=?  1  I 

rVE  g»y  ,  2e»j, 

nVE^i-y^/E^y\  J-Ve  ^vsJiVy^iE^y 

TTl-l  _ 


=  n 


=  n 


fVE  y2 

_  1 


_ 1 _ ^ 

J-VE  ^VEyk^^^y 

\  ^  -I  y  -  » 

l[arc.m(l)-are«n(-l)]]” 


(3.44) 
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Consider  the  last  term  of  the  above  equation, 

!*y/E 


I-Ve 

=  nVEji-yyE^y~^^^  V^-vVEdy 

=  ^arcsin{y/y/E) ^  -  y^/Edy 

Ve 


=  E-2Sf 


- dy 

-^sfEyJX-y^lE  ^ 


So, 


/Ve 
Ve 


r 

Ve  Tr\^y/1  -  y'^/E  ^  2 ' 


At  last,  the  variance  of  ^  is 


If  the  variable  ^  is  a  cosine  function  of  the  phase  0,  say 


(3.45) 


(3.46) 


7/2  ry  P 


y  =  '/Ecos{<f)), 


(3.48) 


and  let  the  variable  ^  be 

n  n 

^  =  Ylyi"^Yl  ^cos{<f>i),  (3.49) 

i=l  t=l 

where  the  P.D.F.  of  </»t  is  /$(^),  the  £?[^]  and  uar[’i']  can  be  computed  in  the 
similarly  manner  as 


£;[^]  =  0, 


(3.50) 
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rn  nE 

mr[^]  =  — .  (3  51) 

For  the  noiseless  asynchronous  CFSK  channel,  the  signal  on  frequency  j  is 
given  by  (3.3).  Thus,  we  have 

var[xj]  =  E[x]\  = 
var[yj]  =  E[yJ\  = 

So,  the  expectation  of  pj  is 

E\Pi]  =  E[x]  +  y]]  =  E\x]]  +  E[y]]  =  p^E.  (3.53) 

Assume  that  there  is  some  phase  </>,  .  And,  let  the  effective  Gaussian  noise  be 
ne//,x  and  with  double-sided  P.S.D.  According  to  large  number 
theorem,  the  signal  on  the  noiseless  asynchronous  channel  can  be  approxi¬ 
mated  by  a  signal  on  the  synchronous  channel  with  Gaussian  noise  when 
IS  large.  And,  the  effective  signal  is  approximated  by 

Xj  =  y/J[;Ecos{(j>j)+n,, 

Vj  =  y/ PjEsin{(j)j)  +  Uy. 

3.7  Appendix  II:Proof  of  Theorem  3 

From  Figure  3.7,  it  is  obvious  that  an  upper  bound  of  the  sum  capacity  of 

the  noiseless  asynchronous  CFSK  channel  without  intensity  information  is 
given  by  I{U\X)  where 


I{U\X)  =  H{U)  -  H{U\X). 


(3.55) 
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Since  m  is  the  number  of  active  frequencies  at  the  output  U,  and  the  proba¬ 
bility  of  m  is  then 

=  (3-56) 


Proof  of  (3.56),  (3.30)  and  (3.31)  can  be  found  in  [21].  And,  (3.26)  is  obvious. 
To  obtain  we  need  to  consider  ,  which  denotes  the  number  of  ways 
that  put  k  balls  into  n  bins  so  that  each  bin  contains  at  least  2  balls.  By 
subtracting  the  cases  that  there  are  bins  that  contain  exactly  2  balls, 
can  be  obtained.  For  the  n  bins  and  k  balls,  if  there  are  j  bins  that  contain 
exact  2  balls,  there  are  (p  nj=i  ways  to  put  the  2j  balls  in  to  the  j 

bins.  For  each  of  these  settlements,  there  are  ^_2j-  ways  to  put  the  others 
balls  into  the  remained  bins.  So,  for  the  cases,  the  number  of  cases  that 
there  are  j  bins  contain  exactly  2  balls  is  given  by 


As  a  result,  is  computed  by 

aS=«S-E<L. 


(3.57) 


(3.68) 


By  substituting  (3.57)  in  (3.58),  (3.32)  is  obtained. 


The  probability  PD\i{d\i)  is  equal  to  the  probability  that  there  are  d  bins  that 
contain  exactly  2  balls  by  putting  T  balls  into  i  bins  without  empty.  The 
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number  of  ways  that  put  T  balls  into  i  bins  without  empty  is  For  the 
cases  that  there  are  s  bins  that  contain  exactly  1  ball,  and  d  bins  contain 
exactly  2  balls,  the  number  of  such  settlements  is 


s!7Vf  ^  T 
t— s,T— s,a 


Then,  PD\i{d\i)  is  given  by 


(3.59) 


PD\iid\i)  = 


a 


(1) 

i,T 


(3.60) 


when  d  <1.  For  d  >  i,  pD\i{d\i)  is  0.  By  substituting  (3.59)  into  (3.60), 
(3.27)  is  obtained. 


Consider  the  model  given  in  Figure  3.7,  the  errors  happen  only  when  the 
double  occupied  frequencies  jump  to  inactive.  Such  error  is  with  the  prob¬ 
ability  q  =  p(0|2),  where  p(0|2)  is  given  in  (3.23).  For  the  reason  that  no 
inactive  frequency  that  can  be  activated,  if  there  are  m  active  frequencies 
at  the  output  U,  then  we  must  have  m  <  i  and  there  are  i  -  m  frequencies 
are  deactivated,  which  are  double  occupied  at  input  X.  So,  there  are  (  j 
possible  cases  that  i  -  m  of  d  double  occupied  frequencies  are  deactivated. 
For  each  of  them,  the  probability  is  g«-(l  -  q)d-ii-m)_  So,  pM|D,/(m|d,z)  is 
given  by  (3.28). 

A  known  fact  before  computing  the  value  H{U\X)  is  that  the  only  error  at 
U  is  caused  by  the  deactivation  of  the  double  occupied  frequencies.  Hence, 
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each  of  the  d  double  occupied  frequencies  is  considered  as  a  separated  binary 
channel  with  transition  probability  q.  Let  Eoid)  be  the  expectation  of  d, 
then  we  have 

H{U\X)  =  h{q)ED{d).  (3.61) 

The  expectation  ED{d)  is  given  by 

T  T 

Enid)  =  ^  o(^  P,(j)Pc|/(<i|i).  (3.62) 

d=0  i=l 

By  substituting  (3.62),  (3.61)  and  (3.56)  into  (3.55),  the  upper  bound  (3.25) 


is  obtained. 


Chapter  4 

Multiuser  Detection  on  the 
CFSK  Channel 

4.1  Introduction 

For  a  multiple  access  system,  the  multiuser  detection  operates  on  the  channel 
output  that  includes  information  from  all  the  users.  For  a  particular  user, 
the  receiver  does  not  treat  other  users  as  unwanted  noise.  On  the  other  hand, 
when  the  receiver  is  making  decision  on  one  of  the  users,  it  takes  information 
from  other  users  as  useful  information. 

As  introduced  in  Section  2.4,  there  are  two  categories  of  multiuser  detec¬ 
tion  for  the  CFSK  system.  The  optimal  approach  takes  all  users  as  a  single 
entity  by  constructing  a  super-trellis,  and  does  MLSD  on  the  super-trellis. 
Due  to  complexity  reasons,  the  optimal  approach  is  not  practical  for  large 
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user  number  applications.  Thus,  suboptimal  approaches  are  required  so  that 
lower  decoding  complexity  can  be  achieved  by  sacrificing  some  system  perfor¬ 
mance.  Iterative  detection  can  provide  suboptimal  performance  by  iterative 
decoding  structure  based  on  each  user’s  component  decoder.  Complexity  of 
the  iterative  detection  is  decided  by  complexity  of  the  component  decoder. 
Two  known  iterative  detections  are  interference  cancellation  and  consensus 
decoder.  Known  results  show  that  there  is  a  large  distance  between  the  per¬ 
formance  of  interference  cancellation  and  that  of  the  optimal  approach  [14]. 
In  this  context,  the  consensus  decoder  will  be  discussed  in  detail. 

Basic  structure  of  the  consensus  decoder  has  been  given  in  Figure  2.8.  Each 
user  has  a  component  decoder  in  the  consensus  decoder.  The  component 
decoder  detects  the  user’s  code  sequence  based  on  the  user’s  code  structure. 
However,  it  does  not  have  any  knowledge  about  other  user’s  code  structures. 
This  is  different  from  the  MLSD  approach,  where  all  users’  code  structmres 
are  considered  together.  Thus,  it  seems  that  the  consensus  decoder  has  no 
difference  from  the  single-user  detection.  In  Figure  2.8,  each  component 
decoder  takes  the  a  priori  information  on  the  symbols  from  other  users  as 
parameters.  These  parameters  are  used  in  computing  the  channel  metric. 
In  this  sense,  the  consensus  decoder  does  multiuser  detection  by  exchanging 
the  a  priori  information  among  users  and  using  the  users’  information  in 
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computing  the  channel  metric. 

In  this  chapter,  the  optimal  multiuser  detection  approach  using  MLSD  is 
introduced  first.  Then,  the  iterative  detection  using  consensus  decoder  is 
discussed.  Based  on  different  codes  and  decoding  algorithms,  the  consen¬ 
sus  decoder  is  classified  into  the  MAP-consensus  decoder  and  the  Turbo¬ 
consensus  decoder.  The  channel  metric  in  the  MAP  algorithm  is  derived, 
which  introduced  high  complexity  to  the  consensus  decoder  when  the  system 

load  is  high.  Thus,  some  simplification  approaches  for  the  channel  metric  is 
discussed  in  the  end. 

4.2  Maximal  Likelihood  Sequence  Detection 

4.2.1  Channel  Metric 

The  optimal  multiuser  detection  operates  on  the  received  signal,  which  can 
be  described  by  the  vector  F  given  in  (2.7).  The  decision  is  made  by  choosing 
a  point  in  the  signal  space  which  is  the  closest  in  Euclidean  distance  to  the 
received  signal  expressed  by  F.  An  example  of  the  signal  space  is  given  in 
Figure  2.3.  The  points  in  the  signal  space  can  be  described  by 


—  (/^O,  ^l,...,/iM-l)) 


(4.1) 
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where  nj  describes  the  number  of  transmissions  on  frequency  fj.  In  Chapter 
3,  fpjipj)  in  (3.2)  gives  a  measurement  of  the  distance  on  the  dimension 
between  the  received  signal  and  the  point  in  the  signal  space.  If  all  the  M 
dimensions  are  considered,  then  the  Euclidean  distance  can  be  computed  by 

£>(•?■./!)  =(^)  He''"  >■  (4.2) 

Thus,  the  channel  metric  for  MLSD  is  given  by  D{F,  p). 

4.2.2  MLSD  on  the  Super- Trellis 

The  channel  metric  given  in  (4.2)  implies  that  the  detection  algorithm  must 
consider  all  users’  code  structures  together.  This  means  that  a  super-trellis 
based  on  all  users’  trellises  should  be  constructed.  In  [14],  method  for  con¬ 
structing  the  super-trellis  code  is  described.  The  input  bits  of  the  super  code 
is  given  by  aligning  all  users  input  bits  together.  State  of  the  users’  trellis  are 
also  aligned  together.  For  the  CFSK  system,  the  output  of  the  super  code 
can  be  expressed  by  the  vector  p  in  (2.2.3),  which  denotes  the  points  in  the 
signal  space. 

For  convenience  of  discussion,  the  example  of  the  2-user  super-trellis  in  Figure 
2.4  is  described  here  again  in  Figure  4.1.  Each  user  uses  a  rate  1/2,  memory 
1  code.  Thus,  the  super-code  has  2  input  bits  and  2  memories.  For  the 
super  trellis,  the  state  number  is  4.  The  branch  output  is  constructed  by  the 
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I 


Figure  4.1;  A  two-user  super-trellis. 


output  of  each  user’s  code.  Since  there  are  4  frequencies,  the  output  symbol 
is  a  vector  of  4  entries.  Each  entry  denotes  the  number  of  users  that  transmit 
on  the  frequency.  Using  the  metric  given  in  (4.2),  a  MLSD  procedure  can 
work  on  the  super-trellis  and  give  the  optimal  detection  results.  Note  that 
MLSD  can  be  implemented  by  the  Viterbi  algorithm. 

A  limitation  of  the  MLSD  approach  is  that  the  complexity  of  the  super-trellis 
IS  extremely  high  when  there  are  a  large  number  of  users.  Though  the  number 
of  input  bits  k  and  memory  m  of  the  super  code  increases  linearly  with  the 
user  number,  the  input  alphabet  and  possible  states  increases  exponentially 
with  k  and  m.  Thus, -complexity  of  the  super-trellis  grows  exponentially 
with  the  number  of  users.  This  is  reason  that  the  suboptimal  detection  using 
iterative  algorithm  is  applied. 
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4.3  Iterative  Detection 

4.3.1  Channel  Metric  for  the  Consensus  Decoder 

In  Figure  2.8,  the  component  decoder  takes  the  signal  F  as  input,  and  makes 
decision  by  choosing  a  point  in  the  signal  space  that  is  the  closest  in  Eu¬ 
clidean  distance  with  F.  A  difference  from  the  MLSD  is  that  the  component 
decoder  behaves  on  the  single  user’s  code  structure.  In  this  case,  each  user 
has  a  single-user  signal  space,  which  is  different  from  that  described  in  Figme 
2.3.  Since  all  users  use  the  same  M-ary  FSK  modulation,  their  single-user 
signal  space  is  an  M-dimensional  space.  The  point  in  the  space  has  the 
coordinate  (0, . . . ,  0, 1, 0, . . . ,  0),  where  the  1  is  in  the  position  of  the  co¬ 
ordinate.  There  is  no  way  to  compute  directly  the  distance  between  a  point 
in  the  single-user  signal  space  and  the  received  signal  that  is  in  the  multiuser 
signal  space.  However,  the  probability  that  a  point  in  the  single-user  signal 
space  leads  to  the  received  signal  can  be  computed.  In  other  words,  the 
channel  transition  probability  p{F\fj^^)  can  be  computed.  This  probability 
can  act  as  a  metric  of  the  distance  between  the  point  in  the  single-user  signal 
space  and  the  received  signal.  Computation  of  this  metric  is  possible  because 
that  the  component  decoder  takes  other  users’  information  on  the  frequencies 


as  parameters. 
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Denote  the  set  of  users  by 


^  —  {Uo,Ui,. . .  ,Ut-i}.  (4.3) 

Assume  user  Uo  chooses  frequency  /j°).  If  user  C/„  chooses  a  frequency 

with  probability  p(/j;)),  choices  of  the  other  T-1  users  amount  to  a  frequency 
combination 


u  \jjQ  7  yJju+i  ’  •  •  •  w  /• 


Then,  the  probability  for  the  frequency  combination  cq  given  that  user  Uq 
transmits  on  frequency  is 


Pi<^\fj?)  =  (4.5) 

Given  that  user  Uo  transmits  and  the  frequency  combination  cq  is  trans¬ 
mitted  by  the  other  T-1  users,  the  multiplicity  factor  Hj  for  each  frequency 
fj  can  be  obtained.  For  the  synchronous  transmission,  the  distribution  of 
the  received  energy  p,-  for  /,-  is  given  by  (3.2),  when  /j®>  and  co  are  known. 
Thus,  the  probability  of  pj  given  and  cq  is 


(4.6) 


If  the  transmission  of  each  frequency  in  the  set  S  (see  (2.4))  is  independent, 
the  probability  of  receiving  frequency  vector  F  given  and  Cq  is 
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Let  all  the  possible  combinations  cq  be  in  the  set  Cq.  Then,  the  probabil¬ 
ity  of  the  received  frequency  vector  F  given  that  Uq  transmits  can  be 
represented  by 


p(f'i/r)=  E  .  (4.8) 

co€Cb 

After  substituting  the  probabilities  P(co|/j“^)  and  P{F\co,  in  (4.5)  and 
(4.7),  the  channel  transition  probability  of  the  CFSK  channel  with  intensity 


information  for  the  synchronous  transmission  is 

pws’) = (f)  e  e  •  •  ■  e  n  n ''  (4-9) 


^1=0  *2=0  ix-i=0  Lw=i 


where  is  the  index  of  the  frequency  for  the  user. 


In  Chapter  2,  the  distribution  of  the  received  energy  pj  for  noisy  asynchronous 
channel  is  not  given.  Assume  the  asynchronous  channel  is  noiseless,  then  the 


If  the  AWGN  is  introduced  to  the  channel  with  double-sided  P.S.D.  No,  then 


the  over  variance  of  pj  can  be  approximated  by  No  -I-  pjE.  When  the  user 
number  is  smaller  than  the  frequency  number,  the  average  A^alue  of  pj  is 
lower  than  1.  Especially,  when  the  spreading  factor  is  high,  average  value  of 


Pj  tends  to  0.  In  this  case,  the  probability  (4.10)  can  be  approximated  by 


pipiKfj?) = 


- TTT - 


(4.11) 
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As  a  result,  the  channel  transition  probability  of  the  CFSK  channel  with 
intensity  information  for  the  asynchronous  transmission  is 


*1=0  12=0  IT- 1=0  Lii=i 


j=0 


4.12) 


4.3.2  The  MAP  Consensus  Decoder 


In  the  transmitters  of  the  CFSK  system,  if  the  convolutional  codes  are  used 
for  the  users,  it  is  possible  for  the  consensus  decoder  to  use  the  MAP  [33] 
algorithm  for  the  component  decoders.  An  important  reason  for  using  the 
MAP  algorithm  is  that  the  component  decoders  need  to  output  the  a  posteri¬ 
ori  probabilities  for  the  symbols.  As  discussed  previously,  these  a  posteriori 
probabilities  are  exchanged  to  other  users’  component  decoders  as  the  a  pri¬ 
ori  probability  in  computing  the  channel  metric  given  in  (4.9)  and  (4.12). 
This  consensus  decoder  is  thus  called  MAP-consensus  decoder. 

In  Figure  4.2,  the  structure  of  the  MAP-consensus  decoder  is  described.  For 
the  first  iteration,  all  users  initialize  the  distribution  of  the  a  priori  prob¬ 
abilities  by  the  uniform  distribution.  In  other  words,  p(/j“^)  is  initialized 
to  i  for  all  users  and  all  frequencies.  After  the  iteration,  the  a  posterior 
probabilities  p{fjJ\F)  can  be  obtained.  Instead  of  making  a  decision  on  the 
output  symbols  instantly,  these  a  posterior  probabilities  are  fed  back  to  other 
users  decoders,  and  are  used  as  the  a  priori  probabilities  in  computing  the 
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Figure  4.2:  The  consensus  decoder  using  the  MAP  algorithm. 
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metric.  So,  the  channel  metric  in  the  next  iteration  is  updated  and 

improved. 

To  obtain  the  a  posteriori  probability  p(/j“^|F),  all  operations  in  the  MAP 
algorithm  are  symbol-based.  Thus,  the  final  decision  is  made  on  the  code 
symbols.  This  is  the  reason  that  the  recursive  systematic  codes  (RSC)  are 
used  m  the  system,  since  we  can  obtain  the  information  bits  directly  from 
the  symbols.  Another  reason  for  application  of  the  RSC  is  that  it  is  easy 
to  avoid  the  catastrophic  codes  when  we  randomly  generate  the  RSC  codes. 
[34]  provides  a  symbol-based  MAP  algorithm  for  decoding  the  RSC  codes. 
Let  the  length  of  the  code  be  r .  The  received  signal  from  start  to  the  current 
stage  is  denoted  by  Fl,  where  t  denotes  the  current  stage.  Denote  the  current 
state  on  the  trellis  by  St,  which  takes  the  value  m.  The  trellis  branches  output 
symbols  Xt  with  the  value  s,  which  belongs  to  the  set  S  in  (2.4).  And,  let  the 
current  received  signal  be  Ft.  The  code  symbol  output  by  Then,  the  forward 
and  backward  recursions  are 

ttt(m)  =  p{St  =  m,  F{) 

_  Dy  IZ»e57»(^t,m  ,m)a(_i(m')  (4.13) 

5Zs6sT'»(A>»n\m)at-i(m')  ’ 

Pt{m)  =  p[F^^^\St  =  m) 

ICm  ls{Pt,m\Tn)at{rn)  ' 


(4.14) 
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The  branch  transition  probabilities  are  given  by 

7s(Ft,  m,  m)  =  p{St  =  m,Yt  =  Ft,  Xt  =  slS't-i  =  m') 

=  p{Yt  =  Ft\Xt  =  s,  St-i  =  m,  St  =  rn)q(Xt  =  s\St-i  =  m, 
R{St  =  =  m), 

(4.15) 

We  can  also  have  the  probability 


5t{m,  m)  =  p{St-\  =rn,St  =  m,  F{) 
=  at{m)')s{Ft,rn  ,m)pt{'rn). 


(4.16) 


Thus,  the  a  posteriori  probability  of  the  symbols  are 


p{Xt  =  s|y;  =  ^)  = 


pjXt  =  s,  F{) 
PiFi)  ‘ 


(4.17) 


By  normalization,  the  value  p{F{)  is  cancelled  out  from  the  above  probability. 
Thus,  the  a  posteriori  probabilities  can  be  written  as 


p{Xt  =  s\Fl)  =  Kt{s) 

=  p{Xt  =  s,Fl) 

=  E(m',m)eB(s)  Oit{m)^s{Ft,  rn,m)Pt{m). 


(4.18) 


The  set  B{s)  denotes  the  set  of  branch  transitions  from  state  m  to  m  with 
output  symbol  s.  The  decision  is  made  by  choosing  a  symbol  s  with  the 
maximum  value  for  At(s).  The  a  posteriori  information  p{Xt  =  s\F[)  for  all 
symbols  is  exchanged  to  other  users  as  the  o  priori  information  p{Xt  —  s) 
for  the  symbols. 
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4.3.3  The  Turbo  Consensus  Decoder 

To  improve  the  error  correction  capability  of  the  CFSK  system,  Turbo  codes 
can  be  applied  in  the  system.  Thus,  both  the  transmitter  and  the  consensus- 
decoder  need  to  be  modified.  Since  the  consensus  decoder  does  symbol-by¬ 
symbol  detection  based  on  the  received  frequency  vector,  the  corresponding 
turbo  encoder  should  also  be  symbol  based.  Robertson  introduced  a  band¬ 
width  efficient  coding  scheme  employing  turbo  codes  [35],  which  is  adapted 
by  the  turbo  consensus  decoder  in  this  context. 

In  a  CFSK  system  using  the  Turbo-consensus  decoder,  the  transmitter  should 
also  use  the  Turbo  code  [36],  which  is  also  called  the  parallel  concatenation 
convolutional  code.  The  encoder  structure  of  Robertson’s  scheme  has  some 
differences  from  the  Berrou’s  encoder,  as  given  in  Figure  4.3.  The  interleaver 
IS  symbol  based.  The  puncturing  is  also  symbol-based.  Thus,  the  systematic 
bits  generated  by  the  second  encoder  will  be  sent  out  with  the  codeword  if 
the  corresponding  symbol  is  not  punctured.  For  this  reason,  the  interleaver 
should  interleave  an  even  indexed  symbol  to  an  even  indexed  position,  and 
interleave  an  odd  indexed  symbol  to  an  odd  indexed  position,  so  that  the  two 
codewords  generated  by  the  two  encoders  with  the  same  systematic  bits  will 
not  be  both  sent  out  or  both  punctured.  After  puncturing,  the  codewords 
will  be  mapped  to  the  corresponding  frequencies,  which  is  the  same  as  we 
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have  discussed  for  the  convolutional  encoder. 


Figiire  4.3:  Symbol  based  Turbo  encoder 


In  the  Turbo-consensus  decoder,  each  user’s  Turbo  iterative  decoder  is  shown 
in  Figure  4.4.  The  received  signal  F  is  fed  into  the  two  component  decoders 
alternatively  according  to  the  puncturing  rules  at  the  encoder.  The  com¬ 
ponent  decoders  are  using  MAP  algorithm.  The  received  signal  is  used  to 
compute  the  channel  metric.  If  the  corresponding  symbol  is  punctured,  the 
metric  is  set  to  1  since  nothing  can  be  obtained  from  the  channel.  The  MAP 
algorithm  is  symbol-based  as  the  one  used  for  the  MAP-consensus  decoder. 
Computation  of  at{m),  I3t{m)  and  7s(Ft,  m,m)  is  according  to  (4.13),  (4.14) 
and  (4.15).  Also,  the  a  posteriori  probability  of  a  symbol  is  obtained  by 
(4.18).  Both  component  decoders  output  p{Xt  =  sj^).  They  are  exchanged 
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to  other  users’  decoders  for  computing  the  metric.  Note  that  the  a  posteriori 
probabilities  generated  by  the  outer  decoder  should  be  exchanged  to  outer 
decoder,  and  the  those  generated  by  the  inner  decoder  should  be  exchanged 
to  the  inner  decoder.  The  a  posteriori  probability  generated  by  outer  decoder 
(or  inner  decoder)  will  also  be  feed  to  the  inner  decoder  (or  outer  decoder)  of 
the  same  user  for  computing  the  extrinsic  information  of  the  corresponding 
systematic  bits  (thus  to  compute  il{m  ,m)).  In  computing  the  extrinsic  in¬ 
formation,  the  input  a  priori  information  of  this  component  decoder  should 
be  t^en  from  the  output  a  posteriori  information  of  this  component  de¬ 
coder,  so  that  the  information  for  the  systematic  bits  is  not  reused  in  later 
Iterations  of  the  Turbo  decoding.  At  last,  the  decision  can  be  made  on  the 
output  symbol  by  both  the  inner  component  decoder  or  the  outer  component 
decoder.  In  this  project,  the  decision  is  made  by  inner  component  decoder. 

4.4  Computational  Complexity  Problems  in 
the  Consensus  Decoder 

4.4.1  Computational  Complexity  of  the  Channel  Met¬ 
ric 

In  computing  the  channel  metric  p(F|/j“))  by  (4.8),  the  summation  is  over 
the  set  C  of  all  possible  frequency  combinations.  Since  each  of  the  T  -  1 
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Figure  4.4:  Symbol  based  Turbo  iterative  decoder 
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users  has  M  possible  choices  for  picking  a  frequency  from  the  set  S,  there 
are  possible  frequency  combinations  in  the  set  C.  Thus,  the  computa¬ 
tional  complexity  of  increases  exponentially  with  the  user  number 

T.  The  extremely  high  computational  complexity  in  highly  loaded  systems 
is  impractical  for  realization  of  the  consensus  decoder.  So,  approximation 
of  the  metric  with  a  reduced  complexity  is  critical  for  the  application  of  a 
CFSK  system. 

Ideally,  the  computational  complexity  of  the  approximated  metric  should  in¬ 
crease  linearly  with  the  number  of  users.  By  observing  the  structure  of  the 
consensus  decoder,  it  is  found  that  most  of  the  frequency  combinations  are  of 
small  probabilities.  For  those  combinations  in  which  all  the  T- 1  users  choose 
small  probability  frequencies,  the  probabilities  of  the  combinations  are  also 
very  small.  Omitting  these  small  probability  combinations  will  not  introduce 
serious  errors  into  the  metric.  So,  one  approach  for  simplifying  the  channel 
metric  is  to  consider  only  the  K{T)  most  probable  frequency  com¬ 

binations  in  the  set  C,  where  K{T)  is  with  polynomial  complexity.  Denote 
the  set  of  these  K{T)  combinations  by  C',  (4.8)  becomes 


(4.19) 
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In  searching  the  most  probable  combinations,  if  the  measurement  is  given  by 
pnsmp=(^y)  (4-20) 

^  u=l  j=0 

the  approach  is  called  “Narrow-sense  most  probable  combinations”  (NSMP). 
If  the  measurement  is 

T-l 

PwsMP  =  (4-21) 

W=1 

the  approach  is  called  “Wide-sense  most  probable  combinations”  (WSMP). 
Generally,  it  is  hard  to  find  the  NSMP’s  except  for  the  brute  force  check, 
which  requires  computation  of  all  the  combinations.  There  are  meth¬ 

ods  in  implementing  the  WSMP’s  by  the  polynomial  complexities.  However, 
the  a  priori  information  is  unavailable  in  the  first  iteration.  Thus, 

the  a  priori  information  given  by  the  first  iteration  is  highly  suboptimal. 
And,  by  exchanging  the  imprecise  a  priori  information,  the  users  introduce 
interference  to  each  other.  This  leads  to  poor  performance  for  later  itera¬ 
tions,  especially  when  the  system  load  is  high.  To  fix  this  problem,  NSMP  is 
appUed  in  the  first  iteration,  since  it  is  possible  to  implement  NSMP  in  the 
first  iteration  with  low  complexity. 

In  the  following  discussion,  the  linear  approaches  for  WSMP  is  discussed 
firstly.  Then,  an  NSMP  approach  for  the  first  iteration  is  given,  which  applies 
the  Viterbi  algorithm  in  searching  the  combinations  with  highest  value  for 
(4.20). 


CHAPTER  4.  MULTIUSER  DETECTION  ON  THE  CFSK  CHANNEL74 

4.4.2  The  /("-Most  Probable  Combinations  Approach 
{KMPC) 

In  the  A'MPC  approach,  the  set  of  the  most  probable  frequency  combina¬ 
tions,  C  ,  IS  of  size  K.  In  C',  the  most  probable  combination  considered  is  the 
one  in  which  all  the  T  -  1  users  choose  their  most  probable  frequencies.  For 
other  combinations  considered,  there  are  T  -  2  users  that  will  choose  their 
most  probable  frequencies.  And,  only  one  user  can  choose  a  less  important 
frequency.  To  decide  which  user  can  choose  the  less  important  frequency,  all 
T  - 1  users  less  important  frequencies  are  queued  in  a  list  according  to  their 
values.  In  constructing  the  first  combination  other  than  the  most  probable 
one,  pick  the  first  frequency  from  the  list,  which  belongs  to  one  of  the  T  -  1 
users.  For  this  user,  the  frequency  from  the  list  will  be  chosen  and  each  of  the 
other  T  -  2  users  will  choose  the  most  probable  frequencies.  For  the  second 
combination  other  than  the  most  probable  one,  pick  the  second  frequency 
from  the  list,  which  belongs  to  one  of  the  T  -  1  users.  For  the  other  T-2 
users,  the  most  probable  ones  will  be  chosen.  Repeat  this  procedure  until 
all  the  K  frequency  combinations  are  obtained.  Generally,  the  number  of 
combinations  K  is  independent  of  the  user  number  T,  though  the  value  of 
K  should  be  adjusted  when  T  increases  so  that  the  performance  of  the  sys¬ 
tem  does  not  decrease.  Since  there  are  M  possible  frequencies  for  each  user. 
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the  list  of  the  less  important  frequencies  can  have  at  most  {T  —  1){M  —  1) 
frequencies.  So,  the  maximal  value  of  K  in  the  fr'MPC  approach  is 

Kma^  =  l  +  {T-1){M-1).  (4.22) 

In  this  case,  M  of  each  user’s  frequencies  are  considered  in  constructing  the 
combinations. 

A  disadvantage  for  this  approach  is  that,  when  one  of  the  users  has  frequency 
probabilities  that  are  much  higher  than  those  of  other  users,  the  resulting 
combinations  will  bias  to  information  from  this  user.  Hence,  the  computed 
metric  is  biased  and  will  degrade  the  performance  of  the  consensus  decoder. 

As  an  example,  an  a  priori  probability  distribution  for  the  4-user  case  is  given 
in  Table  4.1.  4  most  probable  combinations  are  considered  in  this  example  {K 
is  4).  To  compute  p(F|/j°^)  for  user  0,  the  frequency  combinations  obtained 
by  using  the  4MPC  are  given  in  Table  4.2.  Prom  these  results,  it  is  found 
that  four  frequencies  of  user  1  are  chosen  while  only  one  frequency  is  chosen 
for  user  2  and  user  3.  This  is  because  user  1  puts  most  of  the  weights  on 
frequency  0,  1,  3,  and  5.  As  a  result,  the  value  of  P(.P|/(°))  will  bias  to  the 
information  towards  user  1.  In  other  words,  user  1  will  interfere  the  detection 


results  of  other  users. 
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4.4.3  The  A:- Most  Probable  Frequencies  Approach  (A:MPF) 

In  the  kUPF  approach,  for  each  of  the  T  -  1  users,  the  k  most  probable 
frequencies  are  considered  in  constructing  the  combination  set  C\  The  most 
probable  combination  considered  is  the  one  in  which  all  the  T—  1  users  choose 
their  most  probable  frequencies.  For  other  combinations  considered,  there 
^  “  2  users  that  will  choose  their  most  probable  frequencies.  For  the  last 
user,  it  is  possible  for  him  to  choose  the  second,  the  third,  and  up  to  the  k^^ 
less  probable  frequency.  So,  in  the  kUPF  approach,  each  of  the  T  -  1  users 
has  fc  -  1  chances  to  choose  a  frequency  other  than  the  most  probable  one. 

And,  the  relationship  between  K  and  T  is  as  follows 

K{T)  =  1  +  (T  -  l){k  -  1).  (4.23) 

Table  4.1:  An  a  priori  probabilities  distribution  for  the  4-user  case. 


PUo) 

Pifi) 

P(f2) 

Pif3) 

PU4) 

P(h) 

P(f6) 

P(f7) 

Uo 

0.2 
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0.1 

0.06 

0.15 

0.13 

0.17 

0,1 

0.1 

Vz 

0 

0.15 

0.3 

0.22 

0.1 

0.08 

0 

0.15 

Table  4.2:  Combinations  considered  by  the  4MPC  and  2MP  for  user  Uq. 
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Note  that,  when  k  is  equal  to  the  number  of  frequencies  M,  the  A;MPF 
approach  is  the  same  as  the  KMPC  approach  where  K  is  equal  to  Kmax  as 
computed  in  (4.22). 

The  A:MPF  approach  can  keep  fairness  among  all  users  in  computing  the 
metric,  since  k  frequencies  for  each  user  are  considered.  So,  interference 
from  a  particular  user  will  not  happen  for  the  A;MPF  approach.  However, 
the  combinations  obtained  in  this  way  are  not  necessarily  the  most  probable 
ones  over  all  the  possible  combinations  in  the  set  C.  With  the  same  number 
of  combinations  considered,  the  fcMPF  approach  is  not  supposed  to  have  a 
better  performance  then  the  /CMPC  approach. 

When  computing  ^(^1/)°^)  for  user  0  from  the  distribution  given  in  Table 

4.1,  the  frequency  combinations  obtained  by  2MPF  {k  =  2)  is  given  in  Table 

4.2.  It  can  be  observed  that,  though  user  1  puts  most  of  the  weights  on 
frequency  0,  1,  3  and  5,  the  results  are  fair  for  all  the  other  users.  Each  user 
has  considered  two  frequencies  in  this  case. 

4.4.4  The  Narrow  Sense  Most  Probable  (NSMP)  Ap¬ 
proach 

In  the  first  iteration,  no  o  priori  information  is  available.  Usually,  uniform 
distribution  is  used  for  initializing  the  probabilities  pif-f).  Thus,  in  search- 
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ing  the  most  probable  combinations  by  the  measurement  in  (4.20),  the  value 
of  the  first  product  171=^  does  not  make  sense.  To  get  NSMP’s  in  the 

first  iteration  only  needs  to  find  the  combinations  by  the  second  product 


Pdist  — 


-  («)' n .( 

j=0 


No  ) 


(4.24) 


which  measures  the  distance  between  the  combination  and  the  received  sig¬ 
nal.  Define  the  pattern  of  a  frequency  combination  by  the  vector  fi  given  in 
(4.1).  If  the  pattern  for  combinations  with  the  highest  probability  product 
PRice  is  known,  the  number  of  these  combinations  can  be  computed  by 


\  Po  J\  lix 

=  (^-1)!  (r-l-//o)! 

Po^T  —  1  —  /io)!  Px^{T  —  1  —  /xo  —  /Lti)! 

(^  ~  1  ~  -  /^i  -  •  •  •  -  A*a/-2)! 

PM—X^'iT  ~  \  —  /iji/_j)! 

=  (^-1)! 


where  (T  -  1  -  -  . . .  -  /tm-i)!  =  0!  =  1.  If  there  is  some  fij  equals 

T  -  1  for  0  <  j  <  M  -  1,  then  /z*.  will  be  zero  for  k  j.  So,  iV  is  1  in  this 
case.  When  the  number  of  frequency  M  is  greater  than  the  number  of  users 
T,  then  it  is  possible  that  each  user  chooses  a  different  frequency.  In  this 
case,  is  1  for  all  0  <  j  <  M  -  1,  and  N  is  equal  to  (T  -  1)!.  Generally, 
iV  has  a  value  between  1  and  (T  -  1)!.  When  the  spreading  factor  is  of  a 
moderate  value  (between  1  and  10),  N  can  be  small.  In  order  to  explore  the 
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advantage  of  CFSK  in  the  aspect  of  channel  capacity,  practical  applications 
are  interested  in  spreading  factors  between  1  and  10.  At  last,  the  channel 
transition  probability  can  be  computed  as 

p(^|y^(0))  _  NPnsMP  =  FlPwSMpPdisU  (4.26) 

where  Pwsmp  is  the  same  for  all  combinations  in  the  first  iteration. 

In  order  to  find  the  pattern  of  the  combinations  with  the  highest  Pdisu  a  pro¬ 
cedure  that  the  T  —  1  users  make  choice  from  the  M  frequencies  is  simulated. 
Suppose  that  there  are  M  stages.  At  the  stage,  there  are  ^lj  users  choose 
the  frequency  fj.  On  each  stage,  if  there  are  i  users  have  made  their  choices, 
then  the  procedure  is  at  the  state  Si.  Thus,  the  procedure  will  always  start 
from  state  So.  After  each  stage,  there  are  at  most  T  —  1  users  choose  the 
frequency.  So,  state  So  can  connect  to  state  So,  state  Si  and  up  to  state 
St-\.  Similarly,  for  state  Sj,  there  are  at  most  T  —1  —  i  users  choose  the 
frequency.  Hence,  state  Sj  can  connect  to  state  Sj,  state  Sj-j-i  and  up  to  state 
St-1.  For  state  Sj-i,  since  all  the  T  —  1  users  have  made  their  choices,  no 
user  can  choose  the  frequency  in  the  current  stage.  So,  state  Sr-i  can  only 
connect  to  state  St-i. 

With  the  states  and  stages  defined  above,  a  trellis  can  be  constructed  for 
describing  the  frequency  choosing  procedure.  In  Figure  4.5,  a  trellis  of  8- 
user  S-firequency  system  is  described.  The  purpose  is  to  find  the  combination 
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Figure  4.5:  A  trellis  for  searching  the  pattern  of  combinations  with  the  high¬ 
est  Pfiice  in  a  8-user  8-frequency  system. 


pattern  for  computing  (4.26).  So,  we  need  to  consider  the  choice  of  t/j,  U^, 
...,U7.  Thus,  there  are  8  states  and  8  stages  in  Figure  4.5.  At  the  stage, 
if  the  state  is  Si  and  users  choose  frequency  fj,  then  the  next  state  should 
be  The  trellis  start  with  state  since  no  user  make  choose  at  first. 

And,  the  trellis  ends  up  with  state  St-i  since  all  the  T  - 1  users  are  supposed 
to  have  made  their  choices  after  considering  all  the  M  frequencies. 

By  observing  the  trellis  in  Figure  4.5,  it  is  not  diflficult  to  find  that  each 
combination  pattern  denoted  by  p,  can  be  described  by  a  path  through  the 
trellis.  The  purpose  here  is  to  find  a  path  with  the  highest  value  for  Pdist- 
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Pdist  can  be  written  in  another  way  as 


disi 


No 


No 


(4.27) 


Thus,  there  is  a  metric  corresponding  to  each  branch  of  the  trellis.  When 
user  0  does  not  send  frequency  fj,  the  metric  for  the  branch  fij  =  r  on  the 
stage  is  computed  as 


m(j,T)  = 


(4.28) 


If  user  0  send  frequency  /j,  the  metric  for  the  branch  y,j  =  r  on  the  stage 
is  then  computed  as 


m(i,  r)  = 


{  J 


(4.29) 


Since  the  metric  is  defined  for  each  branch  of  the  trellis,  it  is  routine  to  find 
a  path  though  the  trellis  with  the  highest  product  P^ist  in  (4.27)  by  Viterbi 
algorithm. 

Figure  4.6  presents  an  example  of  searching  procedure  on  a  trellis  of  a  4-user 
8-frequency  system.  This  procedme  finds  a  path  through  the  trellis  with  the 
highest  Pdist  so  that  the  channel  transition  probability  P{F\f^^^)  given  in 
(4.26)  can  be  computed. 


At  every  stage  of  the  trellis,  the  branch  metric  is  computed  by  (4.29).  The 
partial  path  product  through  each  input  branch  is  computed  at  each  state. 
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After  comparison  at  every  state,  a  branch  leading  to  the  highest  partial  path 
product  is  survived,  and  all  the  other  input  branches  of  the  state  will  be 
deleted.  Note  that  state  So  has  only  one  input  branch  at  each  stage,  this 
branch  becomes  the  survived  branch  automatically.  Similarly,  after  the  first 
stage,  each  state  has  only  one  input  branch.  It  is  kept  naturally.  On  the 
last  stage,  the  partial  path  products  axe  computed  only  at  state  S^,  and  only 
one  branch  is  survived.  Then,  a  trace  back  procedure  along  the  survived 
branch  is  started  from  the  state  S3  on  the  last  stage,  until  the  state  Sq  on 
the  first  stage  is  reached.  The  path  obtained  by  the  trace  back  procedure  is 
the  one  with  the  highest  product  Pdiaf  Moreover,  the  /i^’s  corresponding  to 
each  branch  of  the  path  compose  to  a  combination  pattern.  In  Figure  4.6, 
the  obtained  pattern  is 

/i  =  {0,1, 1,0, 0,0, 1,0}.  (4.30) 

By  (4.25),  the  number  of  combinations  that  have  the  pattern  in  (4.30)  is  6. 
In  the  first  iteration  of  the  MAP-consensus  decoder,  the  a  priori  probability 
PifV)  is  initialized  by  1/M,  which  equals  to  0.125  in  this  example.  Thus, 
PwsMP  and  P(F|/(o))  can  be  obtained  by  (4.21)  and  (4.26)  respectively. 

For  the  Viterbi  algorithm  operating  on  the  trellis  given  Figure  4.5  and  4.6, 
i+1  operations  axe  needed  at  the  state  5^.  So,  Zti  i  operations  are  needed 
at  each  stage.  Since  there  axe  M  stages,  the  computation  needed  for  the 
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NoUlions: 

I.  Assume  user  Q  tent  frequency  7. 

2  ’ •  denotes  the  muttipitcity  factor  u  is  0: 

*  -  *  denotes  the  muidpltcity  factor  u  is  1; 

* —  *  denotes  the  muttiplicily  factor  u  is  2; 

“ —  ■  denotes  the  muldplicily  factor  u  is  3. 

3.  '  V  *  denotes  the  the  deleted  branches. 

4  The  path  with  arrows  denotes  the  path  with  the  highest  probability. 

Figure  4.6:  An  example  for  the  Viterbi  algorithm  in  a  4-user  8-frequency 
system. 


Viterbi  algorithm  is 


M 


T{T+1) 


(4.31) 


Hence,  the  computational  complexity  of  Viterbi  algorithm  is  0{T^),  which 
is  a  pol)momial  complexity.  This  is  also  the  computational  complexity  for 
NSMP  in  the  first  iteration. 


Chapter  5 

Performance  of  the  Iterative 
Detector 


5.1  Introduction 

Performances  of  the  synchronous  and  asynchronous  systems  are  constrained 
by  the  capacity  results  given  in  Chapter  2.  When  the  spreading  factor  is 
of  moderate  and  high  values,  the  performances  of  the  synchronous  and  the 
asynchronous  system  are  similarly.  In  this  Chapter,  performance  of  the  syn¬ 
chronous  CFSK  system  will  be  discussed  basically.  In  discussing  the  perfor¬ 
mance  of  the  simplified  metrics,  performance  of  the  asynchronous  systems 
will  also  be  given,  since  we  are  more  interested  in  the  applications  of  the 
asynchronous  system  when  the  spreading  factor  is  of  moderate  values.  When 
the  spreading  factor  is  very  low,  i.e.  the  system  load  is  heavy  loaded,  results 
m  Chapter  2  display  that  the  capacity  of  the  asynchronous  system  is  much 


84 


CHAPTER  5.  PERFORMANCE  OF  THE  ITERATIVE  DETECTOR  85 


worse.  Especially,  when  the  spreading  factor  tends  to  0,  the  capacity  of  the 
asynchronous  system  also  approaches  0.  Thus,  we  are  more  interested  in 
application  of  the  synchronous  system  in  this  region.  Results  show  that  the 
current  system  can  not  achieve  decent  performance  in  this  region  even  for  the 
synchronous  system,  since  the  error  floor  is  extremely  high.  Specially  code 
design  schemes  should  be  applied  for  the  system.  This  will  be  the  subject  of 
Chapter  6. 

MLSD  gives  the  optimal  detection  results  for  the  CFSK  system.  Performance 
given  by  iterative  detection  (implemented  by  MAP-consensus  decoder)  is 
suboptimal.  Error  floor  may  exist  for  both  MLSD  and  iterative  detection. 
When  the  codes  are  not  uniquely  decodable,  there  is  error  floor  for  MLSD 
and  iterative  detection.  This  part  of  error  floor  can  be  very  low.  For  itera¬ 
tive  detection^  when  system  load  is  high,  high  error  floor  exists  due  to  the 
suboptimality  of  the  iterative  detection. 

When  the  simplification  approaches,  such  as  AMPC  and  fcMPF,  are  applied 
in  the  iterative  detection,  the  error  floor  is  much  higher.  This  is  because 
that  the  performance  of  AMPC  and  fcMPF  in  the  first  iteration  is  poor. 
Using  NSMP  in  the  first  iteration  can  improve  the  performance  of  simplified 
iterative  detection. 


The  suboptimality  of  the  iterative  detection  is  discussed  in  detail.  Also,  vi- 
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bration  for  the  convergence  of  the  iterations  is  analyzed.  Possible  approaches 
for  solve  the  error  floor  problem  is  discussed.  As  the  first  attempt  to  solve 
this  problem,  channel  interleavers  are  applied  in  the  system.  The  system 
performance  is  improved  when  the  system  load  is  not  high. 

5.2  Comparison  of  MLSD  and  Iterative  De- 
tection 

As  introduced  in  Section  2.4.2,  MLSD  gives  the  optimal  detection  results  for 
the  CFSK  system.  Due  to  the  complexity  of  MLSD,  iterative  detection  is 
preferred,  which  is  a  suboptimal  approach.  For  application  when  system  load 

15  low,  we  expect  that  the  iterative  detection  can  approach  the  performance  of 
MLSD.  When,  system  load  is  higher,  the  suboptimality  of  iterative  detection 
makes  the  performance  much  worse  than  MLSD.  To  observe  the  performance 
of  MLSD,  only  a  couple  of  simple  systems  can  be  simulated  so  that  the 
complexity  of  MLSD  is  not  forbiddingly  high.  Thus,  we  choose  the  codes 
given  in  Figure  2.4  and  2.5  for  simulation.  Each  code  block  has  250  symbols. 

16  iterations  are  used  in  the  iterative  detection. 

In  Figure  5.1,  the  average  bit  error  rate  (BER)  performance  for  the  2-user 
synchronous  CFSK  system  is  shown.  For  each  user,  a  rate  1/2  convolutional 
code  IS  randomly  chosen.  The  number  of  frequencies,  M,  is  2^  =  4.  The 


BER 
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Figiire  5.1:  The  average  BER  performance  of  the  2-user  system  with  codes 
given  in  Figure  2.4. 
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code  set  is  given  in  Figure  2.4.  Since  the  code  set  is  uniquely  decodable,  no 
error  floor  is  observed  in  Figure  5.1  for  the  BER  performance  of  MLSD.  The 
iterative  detection  (ID)  also  provides  a  decent  performance,  which  approxi¬ 
mates  to  the  optimal  performance  given  by  the  MLSD.  The  suboptimality  of 
iterative  detection  does  not  show  up  in  Figure  5.1,  because  the  system  load 
is  not  high  (spreading  factor  is  2). 
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The  average  BER  performance  for  the  3-user  synchronous  CFSK  system  is 
given  in  Figure  5.2,  where  the  codes  are  described  in  Figure  2.5.  Since  this 
code  set  is  not  uniquely  decodable,  the  BER  curve  of  the  MLSD  approaches 
to  an  error  floor  when  the  SNR  is  high.  The  level  of  this  error  floor  is 
decided  by  the  multiplicity  of  the  paths  through  the  super-trellis  with  the 
same  output  sequence.  For  the  iterative  detection,  a  high  level  error  floor 
is  observed.  As  discussed  before,  the  non-unique  decodablity  of  the  codes 
contributes  to  a  part  of  the  error  floor.  The  level  of  this  part  of  error  floor 
equals  that  of  the  MLSD.  In  this  case,  the  system  load  is  higher  (spreading 
factor  is  1.33),  and  the  suboptimality  of  iterative  detection  makes  the  error 
floor  much  higher. 

5.3  Performance  of  the  Simplified  Metrics 

5.3.1  Applications  for  Moderate  Spreading  Factors 

When  the  system  load  is  moderate,  the  simplified  metrics  using  AMPC  or 
A:MPF  approach  may  perform  worse.  This  is  because  that  no  a  priori  infor¬ 
mation  is  available  in  the  first  iteration.  However,  ATMPC  and  fcMPF  will 
depend  on  the  a  priori  information  to  choose  the  most  probable  frequency 
combinations.  The  imprecision  resulted  in  the  first  iteration  is  exchanged 
among  the  users  in  later  iterations.  When  the  system  load  is  moderate  or 
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high,  the  a  prion  information  in  the  first  iteration  is  much  more  imprecise, 
and  such  imprecision  cannot  be  solved  in  later  iterations.  Thus,  poor  perfor¬ 
mance  is  anticipated  for  i^TlVIPC  and  A:MPF  approaches. 

If  NSMP  is  applied  in  the  first  iteration,  the  performance  of  the  first  iteration 
is  supposed  to  approach  that  of  the  optimal  metric,  i.e.  the  metric  without 
simplification.  Thus,  the  a  priori  information  given  by  the  first  iteration 
is  less  suboptimal.  This  makes  it  easier  for  later  iterations  to  approach  to 
better  performance,  though  KMPC  or  fcMPF  is  used  in  later  iterations. 

To  observe  the  effects  discussed  above,  the  synchronous  systems  using  rate 
1/3  codes  with  memory  6  are  simulated.  The  number  of  frequencies,  M,  is 
8.  The  number  of  users  is  varied  from  1  to  4.  The  code  block  length  is  250 
symbols,  and  16  iterations  are  used.  For  the  simplification  approaches,  2 
combinations  is  used  for  i^MPC  and  kMP  when  user  number  is  2.  Thus,  T 
is  2  and  FT  is  2.  From  (4.23),  A:  is  2  for  fcMPF  approach.  Discussion  in  Section 
4.4.3  shows  that  for  the  2-user  system,  FrMPC  and  A:MPF  are  exactly  the 
same  when  K  =  k.  Thus,  2MPC  and  2MPF  have  the  same  performance.  For 
the  4-user  system,  22  combinations  are  used  for  KMPC  and  A:MPF.  Thus, 
T  is  4  and  K  is  22.  Prom  (4.22),  the  maximum  number  of  combinations 
used  for  FfMPC  and  fcMPF,  Kmax,  is  22.  Also,  from  (4.23),  A;  is  8  for  A:MPF 
approach  when  K  is  22.  Since  K  =  Kma.,  KMPC  and  A:MPF  are  exactly 
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that  same.  Thus,  22MPC  and  8MPF  have  the  same  performance  for  the 
4-user  system. 

The  results  are  given  in  Figure  5.3.  The  single  user  performance  provides  a 
lower  bound  for  other  systems  with  the  same  code  rate  and  same  memory. 
The  performance  of  the  2-user  system  with  the  optimal  metric  is  about  O.bdB 
worse  than  the  single  user  performance.  And,  the  performance  of  the  4- 
user  system  with  the  optimal  metric  is  nearly  2dB  worse  than  the  single 
user  performance.  No  error  floor  can  be  observed  for  the  optimal  metric 
performance.  But,  this  does  not  exclude  the  possibility  that  an  extremely 
low  level  error  floor  exists  since  the  codes  are  randomly  generated,  and  the 
unique  decodability  is  not  guaranteed.  When  2MPC  or  2MPF  is  used  for  the 
2-user  system,  the  performance  is  much  worse.  The  BER  curve  converges  to 
a  high  error  floor.  The  error  floor  is  due  to  the  high  suboptimality  of  the 
AMPC  and  A:MPF  approaches.  When  NSMP  approach  is  used  in  the  first 
iteration,  and  2MPC  or  2MPF  is  used  in  later  iterations,  the  performance 
is  much  better,  which  approximates  the  performance  given  by  the  optimal 
metric.  For  the  4-user  system,  similar  results  can  be  observed  for  22MPC 
and  8MPF,  which  have  much  higher  error  floor.  When  NSMP  is  applied  in 
the  first  iteration,  the  error  floor  is  lowered  to  a  much  lower  level.  However, 
performance  given  by  NSMP  plus  22MPC/8MPF  is  still  more  than  IdB 


CHAPTER  5.  PERFORMANCE  OF  THE  ITERATIVE  DETECTOR  92 

worse  than  the  performance  of  the  optimal  metric.  This  is  because  that 
the  spreading  factor  in  this  case  is  2,  which  is  lower  than  the  2-user  system 
(spreading  factor  is  4).  The  suboptimality  for  later  iteration  using  22MPC 
and  8MPF  is  too  high  when  spreading  factor  is  2. 


Figure  5.3:  The  average  BER  performance  of  the  synchronous  4-user  system 
with  rate  1/3  code,  memory  6. 


Since  we  are  more  interested  in  the  asynchronous  system  in  practical  appli¬ 
cations,  we  also  simulated  the  asynchronous  CFSK  system  with  rate  1/3, 
memory  6  codes.  Results  are  given  in  Figure  5.4.  The  performance  given  by 
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the  optimal  metric  is  worse  than  that  of  the  synchronous  system.  Error  floor 
can  be  observed  for  both  the  2-user  system  and  the  4-user  system.  Especially, 
for  the  4-user  system,  high  error  floor  is  observed  for  the  optimal  metric  per¬ 
formance.  Such  degradation  of  performance  is  decided  by  the  capacity  results 
for  asynchronous  system  given  in  Chapter  3.  However,  we  can  still  have  the 
same  conclusions  for  the  simpliflcation  approaches.  When  2MPC/2MPF  is 
used  for  2-user  system,  or  22MPC/8MPF  is  used  for  the  4-user  system,  the 
error  floor  is  much  higher  than  the  one  given  by  the  optimal  metric.  Also, 
application  of  NSMP  in  the  first  iteration  lowers  the  error  floor. 

In  Figure  5.5,  convergence  of  the  asynchronous  4-user  system  is  displayed. 
For  the  first  iteration,  BER  of  22MPC/8MPF  is  much  higher  than  the 
BER  of  the  optimal  metric.  The  poor  performance  in  the  first  iteration 
for  22MPC/8MPF  leads  to  small  step  size  of  convergence  in  later  iterations. 
Thus,  the  iteration  converges  to  a  much  higher  error  level.  When  NSMP 
is  used  in  the  first  iteration,  the  BER  approximated  the  one  given  by  the 
optimal  metric.  Though  22MPC/8MPF  is  used  in  later  iterations,  the  step 
size  of  convergence  is  larger.  Thus,  the  curve  converges  to  a  lower  level. 
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Figure  5.5:  The  convergence  of  the  asynchronous  4-user  system  with  rate  1/3 
code,  memory  6. 
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5.3.2  Applications  for  High  Spreading  Factors 

To  observe  the  performance  of  iterative  detection  when  spreading  factor  is 
high,  more  frequencies  should  be  used  in  the  system.  Thus,  we  use  rate  1/6 
and  2/6,  memory  4  codes  so  that  the  number  of  frequencies,  M,  is  2®  =  64. 
Number  of  users  is  varied  from  1  to  12.  Thus,  the  spreading  factor  is  varied 
from  5.3  to  64.  FVom  results  given  in  Chapter  3,  the  synchronous  system  and 
the  asynchronous  system  perform  similarly.  This  can  be  observed  in  Figure 
5.6,  where  the  synchronous  4-user  system  is  compared  with  the  corresponding 
asynchronous  system.  Spreading  factor  is  16  in  this  case.  4MPC  is  used 
for  the  simulation.  When  BER  is  higher  than  IQ-^,  the  difference  between 
the  two  systems  is  smaller  than  O.bdB.  When  BER  is  higher  than  10-^  the 
difference  is  around  IdB.  Though  a  higher  error  floor  can  be  observed  for  the 
asynchronous  system,  the  difference  between  the  two  system  will  get  much 
smaller  when  the  spreading  factor  gets  higher.  Since  we  are  more  interested 
m  the  application  of  the  asynchronous  system  when  spreading  factor  is  high, 
later  discussions  will  focus  on  the  asynchronous  system.  For  synchronous 
system,  similar  conclusions  can  be  obtained. 

For  rate  1/6  codes  with  memory  4,  it  is  too  complex  to  construct  the  super- 
trelhs  for  MLSD.  Thus,  only  iterative  detection  is  considered  for  the  high 
spreading  factor  applications.  For  iterative  detection,  the  optimal  metric 
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Figure  5.6:  Comparison  of  synchronous  and  asynchronous  4-user  system  with 
rate  1/6  codes,  memory  4. 
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is  also  too  complicated,  since  the  complexity  is  0(64^“^).  Thus,  it  is  only 
possible  to  use  the  simplified  metrics  for  the  high  spreading  factor  applica¬ 
tions.  This  does  not  affect  later  discussions,  since  the  simplified  metrics  can 
achieve  similar  performance  as  the  optimal  metric  when  spreading  factor  is 
high.  This  can  be  observed  in  Figure  5.7,  where  the  asynchronous  4-user  sys¬ 
tem  with  rate  1/6,  memory  4  codes  are  simulated.  The  results  are  given  by 
4MPC,  2MPF  and  NSMP  plus  4MPC.  It  can  be  observed  that  the  difference 
among  the  three  curves  is  quite  small,  which  is  within  O.bdB.  Though  NSMP 
plus  4MPC  seems  to  achieve  a  lower  error  floor,  the  advantage  is  small.  It  is 
reasonable  to  expect  that  the  performance  given  by  NSMP  plus  4MPC  can 
represent  the  performance  given  by  the  optimal  metric  for  the  high  spread¬ 
ing  factor  application.  Hence,  we  can  conclude  that  4MPC  and  2MPF  can 
achieve  similar  performance  given  by  the  optimal  metric. 

In  Figure  5.8,  performances  given  by  4MPC  for  the  asynchronous  1-,  2-,  4- 
,  8-  and  12-user  system  with  memory  4,  rate  1/6  codes  are  given.  It  can 
be  observed  that  the  12-user  system  has  a  high  error  floor,  which  is  around 
4x10  .  For  the  other  4  systems,  the  error  floors  are  all  lower  than  10“^.  Such 
performance  is  decent  for  the  asynchronous  system  with  4MPC  approach. 

For  the  results  given  in  Figure  5.8,  when  there  are  12  users,  the  spreading 
factor  is  5.3,  and  the  sum  rate  of  the  system  is  0.1875  (bits  per  channel 
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Figure  5.7:  Performance  of  asynchronous  4-user  system  with  rate  1/6  codes, 
memory  4. 
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Figure  5.8:  Performance  of  asynchronous  CFSK  system  with  1,  2,  4,  8  and 
12  users,  using  rate  1/6  codes,  memory  4. 
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use  per  chip).  However,  the  channel  capacity  of  the  synchronous  system 
is  around  0.8  (bits  per  channel  use  per  chip)  when  spreading  factor  is  5.3. 
From  results  given  in  Chapter  3,  the  asynchronous  system  has  similar  channel 
capacity  for  this  spreading  factor.  Thus,  the  achieved  sum  rate  is  only  23% 
of  the  channel  capacity.  To  make  use  of  the  channel  more  efficiently,  higher 
code  rate  should  be  applied.  Thus,  the  systems  with  rate  2/6  and  memory  4 
codes  are  simulated.  The  results  are  given  in  Figure  5.9.  The  performances 
for  these  systems  are  worse.  When  user  number  is  8,  the  error  floor  is  around 
2  X  10““*,  while  no  error  floor  can  be  observed  for  the  system  with  rate  1/6 
codes.  When  user  number  is  2  and  4,  the  performance  is  still  decent.  For  the 
8  user  system,  spreading  factor  is  8  and  the  channel  capacity  is  around  0.5. 
The  achieved  sum  rate  is  0.25  (bits  per  channel  use  per  chip)  for  the  8-user 
system  with  rate  2/6  codes,  which  is  50%  of  the  channel  capacity.  Thus,  the 
application  of  the  channel  is  much  more  efficient.  When  user  number  is  12, 
the  performance  of  this  system  is  much  worse  than  the  result  given  in  Figure 
5.8,  which  is  ignored  in  Figure  5.9. 

5.4  Suboptimality  of  Iterative  Detection 

The  suboptimality  of  iterative  detection  may  lead  to  high  error  floor  when 
spreading  factor  is  high.  If  jK'MPC  or  fcMPF  is  applied,  the  iterative  detec- 
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Figure  5.9:  Performance  of  asynchronous  CFSK  system  with  1,  2  4  and  8 
users,  using  rate  2/6  codes,  memory  4. 
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tion  is  more  suboptimal.  Thus,  the  study  of  the  suboptimahty  of  iterative 
detection  is  significant  in  improving  the  system  performance  for  high  spread¬ 
ing  factor  applications. 

At  first,  the  2-user  system  is  considered,  which  is  the  simplest  case.  In  the 
very  beginning,  no  a  prior  information  is  available.  It  is  highly  probable 
that  the  decisions  made  by  the  two  users’  decoders  are  not  correct.  Thus, 
the  computed  a  posteriori  information  for  the  symbols  is  imprecise.  When 
the  information  is  exchanged  to  the  other  user,  it  is  equivalent  to  introducing 
the  multiuser  interference  (MUI)  to  the  other  user. 

For  the  next  iteration,  whenever  the  component  decoder  cannot  get  over  the 
MUI,  new  errors  are  generated  in  computing  the  a  posteriori  information. 
So,  the  MUI  still  exists  and  is  delivered  to  the  other  user  again. 

If  the  MUI  is  not  resolved  after  many  iterations,  the  error  floor  of  the  iterative 
detection  is  formed.  There  are  four  cases  that  the  MUI  can  influence  the  error 
floor: 

1.  If  both  component  decoders  are  powerful,  the  MUI  can  be  resolved 
by  the  decoders.  Then,  there  is  no  error  floor  caused  by  the  iterative 
detection. 


2.  If  one  of  the  component  decoders  is  powerful,  it  can  resolve  the  MUI 


CHAPTER  5.  PERFORMANCE  OF  THE  ITERATIVE  DETECTOR  104 

whenever  it  gets  the  MUI.  However,  the  other  decoder  cannot  resolve 
the  MUI.  The  error  floor  will  get  lower  and  lower,  and  disappear  finally. 

Such  convergence  is  slow  since  only  one  component  decoder  work  prop¬ 
erly. 

3.  If  one  of  the  component  decoders  is  powerful  and  can  resolve  the  MUI 
it  obtains.  The  other  component  decoder  is  not  powerful,  and  sends 
the  MUI  to  the  powerful  component  decoder.  However,  the  MUI  is  too 
high,  and  the  powerful  component  decoder  is  not  strong  enough  to  get 
rid  of  the  MUI.  Thus,  such  MUI  is  delivered  between  the  two  decoders. 
This  will  lead  to  the  vibration  of  the  convergence.  Such  vibration  can 

be  observed  in  Figure  5.5  for  the  convergence  of  8MPC  after  the  10 
iterations. 

4.  If  both  component  decoders  are  not  powerful,  the  MUI  is  exchanged 

between  the  two  component  decoders.  The  error  floor  is  high  in  this 
case. 

Several  factors  can  affect  the  level  of  the  MUI: 

1.  In  the  first  iteration,  the  decision  is  not  precise.  MUI  are  generated. 
If  correlation  of  the  two  users’  codes  is  strong,  the  MUI  will  be  much 


more  serious. 
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2.  In  later  iterations,  the  correlation  of  the  users’  codes  will  still  aflfect  the 
level  of  the  MUL 

3.  For  all  iterations,  the  error  correcting  ability  of  the  codes  will  affect 
the  level  of  the  MUL  If  the  codes  are  powerful,  the  MUI  are  resolved. 
Otherwise,  the  MUI  are  kept  or  may  be  more  serious. 

4.  If  the  simplified  metrics  are  used,  the  MUI  will  be  more  serious  or  new 
MUI  is  generated. 

For  the  system  with  more  than  two  users,  the  exchange  of  MUI  is  much  more 
complicated.  But,  similar  relationship  can  be  obtained  between  MUI  and 
error  floor. 

Though  the  convolutional  codes  applied  in  the  previous  discussions  can  be 
very  powerful,  when  the  coded  sequence  are  mapped  to  the  M-ary  symbols, 
the  Hamming  distance  of  the  codes  is  decreased.  For  example,  the  4-tuple 
1100  are  mapped  to  symbol  sequence  30  (4  frequencies  used),  and  0011  are 
mapped  to  03.  The  Hamming  distance  of  the  bit  sequence  is  4,  the  Hamming 
distance  of  the  4-ary  symbol  sequence  is  only  2.  Thus,  the  signal  mapping 
in  this  way  may  degrade  the  code  performance  a  lot.  This  is  the  reason  that 
high  error  floor  is  formed  when  there  are  too  many  users. 


To  reduce  the  suboptimality  of  iterative  detection,  there  are  3  possible  ap- 
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proaches: 

1.  Decrease  the  correlation  of  the  codes  by  using  the  channel  interleavers. 
When  there  are  many  users  in  the  system,  the  MUI  is  too  high  after  the 
first  iteration.  Thus,  the  channel  interleavers  can  improve  the  system 
performance  only  when  spreading  factor  is  of  moderate  or  high  values. 

2.  In  order  that  the  MUI  is  not  serious  after  the  first  iteration  if  A'MPC  or 
fcMPF  is  used,  NSMP  approach  should  be  applied  in  the  first  iteration. 

3.  By  using  non-binary  codes,  such  as  the  Reed-Solomon  codes  or  con¬ 
volutional  code  over  GF{q),  the  distance  property  of  the  codes  can  be 
guaranteed.  Thus,  the  codes  are  more  powerful  in  decreasing  the  MUI 
level. 

For  the  above  three  approaches,  NSMP  is  already  discussed  in  Section  4.4.4. 
The  application  of  channel  interleavers  will  be  discussed  in  the  following 
section.  Code  design  issues  will  be  given  in  Chapter  6. 

5.5  Effects  of  Channel  Interleavers 

The  channel  interleaver  is  symbol  based,  which  permutes  each  n-tuple  coded 
bits  randomly.  As  denoted  in  Figure  2.1,  the  channel  interleaver  is  con¬ 
catenated  after  the  encoder  at  the  transmitter.  At  the  consensus  decoder. 
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the  users’  component  decoders  also  need  modifications,  as  denoted  in  Figure 
5.10.  Thus,  each  user’s  component  needs  to  know  the  other  T  —  1  users’ 
interleaver  structures.  This  makes  the  structure  of  the  component  decoder 
more  complicated.  However,  the  computational  complexity  of  the  decoding 
algorithm  does  not  increase. 


Figure  5.10:  The  user’s  component  decoder  with  channel  interleavers. 


The  interleaver  size  can  affect  the  performance  of  the  detection.  Generally, 
the  interleaver  size  equals  the  length  of  the  coding  block.  Figure  5.11  shows 
the  performance  of  the  synchronous  3-user  system  with  codes  given  in  Fig¬ 
ure  2.5,  and  with  different  interleaver  sizes.  These  results  axe  compared  with 
those  given  in  Figure  5.2.  We  can  observe  that  when  interleaver  size  in¬ 
creases,  the  system  performance  increases.  This  is  because  that  the  random 
channel  interleaver  reduces  the  level  of  correlation  among  the  users’  codes. 
When  the  interleaver  size  is  large  enough,  the  users’  code  will  be  independent 
with  each  other.  In  such  condition,  the  increase  of  interleaver  size  will  not 
improve  the  performance  any  more.  This  can  be  observed  in  Figure  5.11  that 
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when  the  interleaver  size  increases  from  1000  to  10000,  the  performance  only 
improves  a  little,  compared  with  the  large  improvement  when  the  interleaver 
size  increases  from  250  to  1000.  Another  effect  is  that  the  performance  of  the 
system  with  channel  interleavers  is  even  better  than  the  performance  given 
by  MLSD  when  the  interleaver  size  is  large  enough.  This  is  not  surprising, 
since  MLSD  is  without  channel  interleavers,  and  the  interleavers  can  increase 
the  memory  of  the  super-trellis.  Thus,  such  comparison  is  not  fair  for  MLSD. 

Since  the  channel  interleavers  can  improve  the  system  performance,  it  may 
be  interesting  to  observe  the  effects  when  channel  interleavers  are  applied 
for  4-user  system  with  rate  1/3  codes.  As  given  in  Figure  5.3  and  5.4,  poor 
performance  is  obtained  when  4MPC  or  2MPP  is  used.  When  NSMP  is  used 
for  the  first  iteration,  the  performance  is  still  worse  than  that  of  the  optimal 
metric.  After  utilizing  the  channel  interleavers  in  this  system,  the  results  are 
compared  with  those  given  in  Figure  5.4.  In  Figure  5.12,  it  can  be  found 
that  the  performances  of  4MPC/2MPF  and  NSMP  plus  4MPC  are  improved 
when  channel  interleavers  are  used.  For  high  SNR’s,  their  performances 
approximate  the  one  given  by  the  optimal  metric  with  channel  interleavers. 
And,  the  error  floor  with  channel  interleavers  is  lower  than  the  error  floor 
without  channel  interleavers.  In  Figure  5.13,  convergence  behavior  of  the 
system  with  channel  interleavers  is  described.  When  channel  interleavers 
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Figure  5.11:  Performance  of  synchronous  3-user  system  with  channel  inter- 
leavers.(the  codes  are  given  in  Figure  2.5.) 
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are  applied  with  4MPC,  the  performance  of  the  first  iteration  is  still  worse, 
since  the  interleavers  do  not  help  provide  a  priori  information  in  the  first 
Iteration.  However,  the  interleavers  make  the  users’  codes  less  correlated 
with  each  other.  This  help  reduce  the  MUI  among  users.  Thus,  it  is  easier 
for  the  component  decoders  to  reduce  the  MUI  in  later  iterations.  In  Figure 
5.13,  we  can  find  that  the  convergence  curve  for  4MPC  with  interleavers 
goes  down  with  a  larger  step  size  than  that  of  the  curve  for  4MPC  without 
interleavers.  Such  effect  can  also  be  observed  for  the  curves  for  the  NSMP 
plus  4MPC  and  the  optimal  metric.  Thus,  the  whole  procedure  converges  to 

a  low  BER  level.  And,  this  make  NSMP  plus  4MPC  performs  more  similar 
to  the  optimal  metric. 

5.6  Summary 

In  this  Chapter,  the  performance  of  iterative  detection  is  discussed,  where 
iterative  detection  is  implemented  by  MAP-consensus  decoder.  For  Turbo¬ 
consensus  decoder,  similar  conclusions  can  be  obtained.  The  iterative  detec¬ 
tion  is  suboptimal  when  compared  with  MLSD,  which  provides  the  optimal 
performance  for  multiuser  detection  of  the  CFSK  channel.  The  suboptimal- 
ity  of  iterative  detection  is  because  that  the  multiuser  interference  (MUI) 
is  exchanged  among  users,  and  the  users’  component  decoders  may  not  be 


Figure  5.12:  Performance  of  asynchronous  4-user  system  with  channel  inter¬ 
leavers,  code  rate  1/3,  memory  6. 
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powerful  enough  to  resolve  the  MUI.  Thus,  error  floor  is  resulted  when  MUI 
cannot  be  reduced  by  the  users’  component  decoders. 

To  lower  the  level  of  the  error  floor,  possible  approaches  include  the  appli¬ 
cation  of  NSMP  in  the  first  iteration,  the  application  of  channel  interleaver 
and  the  application  of  powerful  non-binary  code  for  the  users’  codes.  The 
results  have  shown  that  NSMP  can  improve  the  performance  of  /irMPC  or 
fcMPF  effectively.  However,  such  improvement  is  hmited,  especially  when 
the  system  load  is  high.  Channel  interleavers  can  help  improve  the  perfor¬ 
mance  further.  For  moderate  spreading  factors,  channel  interleavers  make 
the  simplified  metrics  work  as  well  as  the  optimal  metric.  Also,  the  error 
floor  of  the  performance  is  lowered.  However  the  help  provided  by  channel 
interleavers  is  limited,  since  the  ability  to  reduce  the  MUI  is  finally  decided 
by  the  error  correction  ability  of  the  channel  codes.  When  spreading  factor 
is  low,  weak  codes  cannot  resolve  the  MUI  problem  at  all.  In  such  condition, 
the  application  of  channel  interleavers  does  not  make  sense. 

An  important  conclusion  is  that  the  key  factor  of  lowering  the  error  floor  is 
the  code  design.  In  previous  discussion,  we  find  that  the  randomly  gener¬ 
ated  binary  convolutional  codes  do  not  provide  decent  performance  since  the 
signal  mapping  decreases  the  Hamming  distance  of  the  codes.  Thus,  codes 
in  the  GF{q)  should  be  considered.  This  will  be  the  subject  of  the  following 
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Chapter. 


Chapter  6 

Code  Design  for  the  CFSK 
System 

6.1  Introduction 

In  former  discussions,  the  channel  codes  are  randomly  generated  RSC.  For 
the  signal  mapping  scheme,  each  n-tuple  coded  bits  are  mapped  to  a  fre¬ 
quency.  As  a  result,  the  symbol  sequence  can  be  considered  as  a  randomly 
chosen  non-binary  trellis  coded  sequence.  Generally,  the  unique  decodability 
of  the  non-binary  code  is  not  guaranteed.  For  applications  of  low  spreading 
factors,  the  number  of  users  is  high  and  the  problem  of  uniquely  decodable 
code  becomes  much  more  serious.  Code  design  for  the  non-binary  trellis 
codes  with  unique  decodability  is  extremely  difficult.  Known  uniquely  de¬ 
codable  block  codes  (UDBC)  provide  us  solution  for  the  design  of  uniquely 
decodable  non-binary  trellis  codes  [17] [37] [38] [39] [28].  In  [17],  Chang  and 
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Wolf  proposed  the  uniquely  decodable  block  codes  with  sum  rates  approach¬ 
ing  the  capacity  bounds  for  T-user  JW-frequency  channel  with  and  without 
intensity  information.  In  [37],  Kautz  and  Singleton  found  the  KS  code  that 
is  based  on  a  maximum-distance-separable  (MDS)  outer  code  (e.g.,  Reed- 
Solomon  code)  and  an  orthogonal  weight-one  inner  code.  The  KS  code  is 
proved  to  be  uniquely  decodable  on  the  binary  adder  channel.  The  orthogo¬ 
nal  weight-one  inner  code  of  KS  code  is  similar  to  the  signal  mapping  function 
m  the  CFSK  system.  In  [39],  example  of  the  appUcation  of  KS  code  based 
on  Reed-Solomon  code  over  GF(7)  is  discussed.  This  provides  us  a  clue  for 
applying  the  Reed-Soloman  codes  in  CFSK  system  for  signal  mapping.  Note 
that  apphcation  of  the  KS  code  in  CFSK  will  lead  to  low  code  rate,  which 
is  far  from  the  capacity  bounds  given  in  Chapter  3.  A  major  problem  for 
the  above  UDBC  is  that  these  codes  are  weak  in  the  error  correction  capa- 
biUty.  A  possible  approach  to  solve  this  problem  is  proposed  in  [28],  where 
the  UDBC  can  be  concatenated  with  the  binary  convolutional  code  so  that 
decent  performance  and  unique  decodability  is  provided.  In  Section  6.2,  the 

appUcation  of  UDBC  given  in  [17]  are  discussed.  Results  of  the  application 
is  given  in  Section  6.3. 
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6.2  Application  of  Uniquely  Decodable  Block 
Codes(UDBC) 

6.2.1  UDBCs  with  Sum  Rate  Approaching  the  Capac¬ 
ity  Bound 


In  [27],  Chang  and  Weldon  proposed  an  iterative  construction  for  the  UDBC 
on  the  binary  memory  less  T-user  adder  channel.  Chang  and  Wolf  extended 
the  construction  to  non-binary  UDBC  on  the  T-user  M-frequency  channel 
with  intensity  information  in  [17]. 


The  iterative  construction  works  on  the  difference  matrix  Dj  in  the  it¬ 
eration.  The  difference  matrix  is  initialized  by  Dq  =  [1].  Then,  for  the  j** 
iteration,  the  difference  matrix  is  given  by 


Dj-i  Dj^i 
Dj—i  Dj—i  , 

/,_i  0 


(6.1) 


where  Dj-\  is  the  difference  matrix  in  the  iteration,  is  the  identity 

matrix,  and  0  is  the  zero  matrix.  By  the  above  iterative  procedure,  Dj 
has  Nj  columns  {Nj  =  2^)  and  Pj  rows  {Pj  =  {j  -|-  2)2-^“^).  For  a  set  of 
codes  constructed  by  j  iterations,  let  D  =  Dj,  then  the  difference  matrix  for 
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b 

ab 

a?b 

a^-^b 

where  a  is  an  arbitrarily  chosen  integer  that  satisfies 


a>Pj{M-l)  +  l. 


Thus,  there  are 


T={M-l){j  +  2)2’-^ 


rows  in  the  matrix  O*.  Then,  each  row  of  D'  can  construct  a  code  for  a 
particular  user.  Hence,  a  set  of  codes  for  T  users  can  by  constructed.  The 
code  has  y  =  2  codewords,  denoted  by  w,  and  w'.  Denote  the  position 

in  wi  and  by  and  01,  ^.  Then,  w,  ,,  and  tcj  ^  can  be  constructed  from 
the  (i,  position  of  b*,  di^^,  as 

I  Wi,k  =  =  0,  if  di^k  =  a', 

<  Wi,k  =  0,  w'i,^  =  o^  =  (6.5) 

V  Wi^k  —  a»,  =  0*2,  if 

Chang  and  Wolf  had  proved  in  [17]  that  the  above  construction  gives  us  a 
set  of  uniquely  decodable  codes  for  T  users  M  frequencies.  For  the  T  users, 
the  sum  rate  can  be  given  by 


-U  +  2), 


(6.6) 
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where  y  denotes  the  number  of  codewords  in  the  code,  and  j  denotes  the 
number  of  iterations  used.  It  is  not  difficult  to  show  that  when  j  tends  to  in¬ 
finity,  the  sum  rate  Rsum  will  approach  the  capacity  bound  of  the  CFSK  chan¬ 
nel  with  intensity  information.  The  sum  rates  of  some  constructed  UDBC 
are  compared  with  the  capacity  bounds  of  the  CFSK  channel  with  intensity 
information.  The  results  are  given  in  Table  6.1. 


Table  6.1:  Comparison  between  sum  rates  of  Chang  and  Wolf’s  UDBC  and 
capacity  results. 


M 

T 

j 

N 

Rsum 

percent  of  capacity 

6 

1 

2 

3 

69.74  % 

3 

16 

2 

4 

4 

69.97  % 

3 

224 

5 

32 

7 

73.50  % 

3 

12,288 

10 

1,024 

m 

78.42  % 

m 

23,068,672 

Eg 

1,048,576 

84.05  % 

a 

18,432 

10 

1,024 

18 

76.94  % 

4 

34,603,008 

20 

1,048,576 

33 

83.11  % 

For  the  6-user  3-frequency  case,  1  iteration  is  used.  Thus,  D  can  be  con¬ 
structed  from  Do  =  [1]  as 


D  =  Di  = 


1  1 

1  -1 

1  0 


Let  the  integer  a  be 


(6.7) 


a  =  Pi{M  -  1)  +  1  =  7, 


(6.8) 
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where  Pi  equals  3.  Thus,  the  difference  matrix  D*  is 


D*  = 


1 

1 

1 

7 

7 

7 


1 

-1 

0 

7 

-7 

0 


(  Co 
Cl 
C2 

C3 
C4 
Cs 


(6.9) 


Prom  D*,  the  code  set  for  six  users  is  constructed  as 

{(1.1).  (0,0)},  } 

{(1.0),  (0,1)}, 

{(1.0),  (0,0)}, 

{(7, 7),  (0,0)},  f-  (6.10) 

{(7,0),  (0,7)}, 

{(7,0),  (0,0)}  ^ 

To  apply  this  code  set  in  the  CPSK  system,  we  can  use  the  following  map 

(6.11) 

As  a  result,  we  obtain  the  non-binary  UDBC  for  the  CFSK  channel  with 
intensity  information  as 


'  Oo  =  {(/l, 1 

Cl  =  {(/i,/o),(/o,/i)}, 

,  —  {(/l./o),  (/o, /o)}, 

C3  =  {(/2,/2).(/o./o)}, 

C4  =  ((f2,M,(/o,/2)}, 

I  =  {(/2./o),(/o,/o)}  , 


(6.12) 


6.2.2  Concatenation  with  Convolutional  Codes 


It  is  easy  to  prove  that  the  code  set  given  in  (6.12)  does  not  have  decent  error 
correction  performance.  To  endow  the  system  with  strong  error  correction 
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capability,  it  is  necessary  to  use  error  correction  code  with  better  performance 
before  applying  the  UDBC.  Thus,  convolutional  codes  are  utiUzed  in  the 
system  before  signal  mapping  (using  the  UDBC).  The  transmitter  described 
in  Figure  2.1  is  modified  as  described  in  Figure  6.1.  As  discussed  previously, 
the  symbol-based  interleaver  is  optional.  The  coded  sequence  output  by 
the  convolution  code  is  then  encoded  by  the  UDBC.  Since  there  are  two 
codewords  in  each  user’s  code,  each  codeword  represents  1  bit.  Thus,  each 
bit  of  the  coded  sequence  output  by  the  convolutional  code  is  mapped  to 
the  a  codeword  in  the  UDBC  set.  In  Figure  6.1,  an  example  is  given  by  the 
code  Co  for  user  Uq.  The  encoding  rule  for  UDBC  is  that  0  is  mapped  to 
the  codeword  (1, 1)  and  1  is  mapped  to  the  codeword  (0, 0).  For  the  M-FSK 
modulation,  the  S3rmbol  1  is  than  mapped  to  frequency  fi  and  0  is  mapped 
to  frequency  /g. 


u 

Convolutional 

V 

Symbol-based 

vMOI...^ 

UDBC 

c={noo...}^ 

M-FSK 

Code 

Interieaver 

Figure  6.1:  Convolutional  code  concatenated  with  UDBC  in  CFSK  trans¬ 
mitter. 


At  the  receiver  side,  the  user’s  component  decoder  has  the  same  structure 
as  we  have  discussed  in  Chapter  4.  In  the  MAP  algorithm,  the  symbol  is 
defined  by  the  branch  output  on  the  convolutional  code’s  trellis.  In  other 
words,  the  n-bit  symbol  can  be  represented  by  n  codewords  of  the  UDBC. 
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Computation  and  decision  for  the  MAP  algorithm  is  made  on  the  symbols, 
since  we  need  the  a  posteriori  probabilities  of  the  symbols.  Furthermore, 
information  bits  can  be  obtained  from  decisions  of  the  symbols,  since  RSC  is 
used  for  convolutional  code  in  the  transmitter.  Note  that  UDBC  endows  the 
trellis  code  with  unique  decodability,  we  can  use  the  same  RSC  for  all  users. 
For  the  (n,k)  RSC,  each  n-tuple  of  the  coded  bits  is  mapped  to  a  symbol. 
Thus,  there  are  Ng  =  2"  symbols.  Each  symbol  of  user  14  is  denoted  by 

~  /i2,15  •  •  •  >  /jAT,!,  .  .  .  ,  ,  fjN,n}-  (6.13) 

To  implement  the  MAP  algorithm,  the  channel  metric  needs  to  be  modified 
as 


(6.14) 

where  N  denotes  the  length  of  the  UDBC  codeword,  denotes  the  number 

of  transmissions  on  frequency  fj  in  the  position  of  the  UDBC  codeword 


for  the  bit  of  the  symbol.  The  received  signal  F  is  a  M  x  matrix 

Po,N,l  .  •  .  Po,l,n)  P0,2,n)  •  •  •  Po,N,n 

PW  •  •  •  Pl,l,n,  Pl,2,n,  .  .  .  pl,N,n 

PM-1,N,1  •  .  .  PA/-l,l,n)  PM-\,2,ni  •  ■  •  PM-l,N,n 

(6.15) 

m  which  the  element  pj^k,x  denotes  the  received  energy  level  of  frequency  fj 
in  the  fc*'*  position  of  the  UDBC  codeword  for  the  x*'*  bit  of  the  symbol. 


Po,i,i,  P0,2,l> 
P  _  Pl.1,1)  Pi, 2,1, 

.  PM-1,1,1,  PM-1,2,1, 
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Base  on  the  symbol  the  structure  of  the  MAP-consensus  decoder  or  the 
Turbo-consensus  decoder  is  preserved,  and  the  simplification  approaches  for 
the  metric  given  in  Chapter  4  are  also  effective. 

6.2.3  Performance  Results 

To  observe  the  performance  of  the  CFSK  system  using  UDBC,  the  concate¬ 
nation  scheme  given  in  previous  section  is  adopted.  3  frequencies  are  used  in 
simulations.  1  or  2  iterations  are  applied  in  generating  the  UDBC  code  set. 
As  given  in  Table  6.1,  the  maximum  number  of  users  in  these  cases  is  6  or  16. 
Then,  the  spreading  factors  of  applications  for  these  cases  are  generally  less 
then  1.  Prom  results  in  Chapter  3,  the  capacity  of  asynchronous  CFSK  chan¬ 
nel  with  intensity  information  decreases  a  lot  when  spreading  factor  is  less 
than  1,  which  implies  that  only  synchronous  system  is  applicable  for  achiev¬ 
ing  high  coding  rate.  Thus,  synchronous  system  is  used  in  the  simulations. 
For  the  convolution  code,  rate  1/2  RSC  is  used.  Thus,  the  sum  rates  achieved 
by  the  system  is  1/2  of  the  sum  rates  given  in  Table  6.1.  The  length  of  the 
coding  block  is  250  symbols.  Both  channel  interleavers  and  NSMP-fKMPC 
are  applied  for  the  system. 

Results  are  given  in  Figure  6.2.  When  1  iteration  is  used  in  generating  the 
UDBC  codes  (block  length  is  2),  and  6  users  access  to  the  system,  BER  is 


BER 
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Figure  6.2:  Performance  of  3-frequency  system,  using  rate  1/2  convolutional 
codes  concatenated  with  uniquely  decodable  block  codes. 
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below  10  ^  when  SNR  is  higher  than  13dB.  No  error  floor  is  observed  for  the 
performance.  When  2  iterations  are  used  for  UDBC  codes  (block  length  is 
4),  14  users  access  to  the  system,  BER  is  around  10“®  when  SNR  is  13dB. 
Also,  no  error  floor  is  observed.  Such  performances  is  decent  for  practical 
applications.  However,  for  the  2  iteration  case,  maximum  number  of  users 
is  16.  When  user  number  is  larger  than  14,  the  performance  gets  worse.  As 
given  in  Figure  6.2,  the  error  floor  is  extremely  high  (BER  is  around  10“^) 
when  15  users  access  to  the  system.  Also,  it  is  easy  to  predict  that  when  16 
users  access  to  the  system,  high  error  floor  will  be  observed  either. 

Generally,  the  maximum  number  of  users  need  to  access  to  the  system  so  that 
the  maximum  sum  rate  is  achieved.  However,  a  large  amount  of  the  simu¬ 
lation  results  (not  shown  in  Figure  6.2)  display  that  most  of  the  maximmn 
sum  rates  given  in  Table  6.1  are  not  achievable,  where  high  error  floor  can  be 
observed.  The  reason  for  the  high  error  floor  has  been  discussed  in  Chapter  5. 
The  error  correction  capability  of  the  channel  codes  decides  the  level  of  error 
floor  to  a  large  extent,  since  the  channel  codes  are  extremely  important  in 
reducing  the  MUI  among  users.  For  our  simulations,  rate  1/2  RSC  codes  are 
applied,  which  is  not  strong  enough  to  resolve  the  MUI  problem  when  a  large 
number  of  users  access  to  the  system.  Thus,  powerful  error  correction  codes, 
which  is  suitable  for  resolving  the  serious  MUI  in  the  CFSK  system,  are 
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necessary  to  achieve  higher  code  rate.  Nonbinary  codes  such  Reed-Solomon 
codes  or  convolutional  codes  over  GF{q)  [40]  are  candidates  for  substituting 
the  binary  convolutional  codes  in  former  system. 


Chapter  7 

Conclusions  and  Future  Work 

7.1  Conclusions 

Iterative  multiuser  detection  for  CFSK  system  is  investigated  in  this  report. 
The  theoretical  background  is  provided  in  [17]  [18],  which  addressed  the  ca¬ 
pacity  issues  for  the  synchronous  CFSK  channels.  In  Chapter  3,  more  capac¬ 
ity  results  for  the  asynchronous  CFSK  channels  are  discussed.  The  results 
show  that  the  asynchronous  CFSK  channel  perform  similarly  when  the  sys¬ 
tem  load  is  low.  When  the  system  load  gets  higher,  performance  of  the 
asynchronous  CFSK  channels  gets  worse.  Thus,  asynchronous  CFSK  system 
is  suitable  for  application  of  moderate  or  low  system  load  applications.  For 
high  system  load  applications,  only  synchronous  CFSK  system  is  usable.  Re¬ 
sults  given  in  [18]  show  that  the  capacity  of  CFSK  system  is  higher  than  the 
CDMA  system.  When  the  spreading  factor  is  of  moderate  or  high  values. 
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the  advantage  of  capacity  for  CFSK  is  the  greatest.  This  is  the  reason  that 
most  of  our  research  is  focus  on  the  applications  for  both  synchronous  and 
asynchronous  CFSK  system  when  spreading  factor  is  moderate.  We  are  also 
interested  in  the  application  of  synchronous  system  when  spreading  factor 
is  lower  than  1.  However,  it  is  much  more  difficult  in  this  region  to  achieve 
higher  code  rates  than  the  bounds  given  by  the  multiuser  CDMA  system. 

Techniques  for  multiuser  detection  is  discussed  in  Chapter  4.  The  optimal 
detection  results  is  provided  by  the  MLSD  approach.  However,  the  high 
complexity  of  MLSD  is  impractical  for  applications.  Iterative  detection  ap¬ 
proaches  are  preferred.  The  multiuser  channel  metric  given  in  Chapter  4  is 
of  non-exponential  complexity,  which  is  extremely  high  when  the  number  of 
users  or  frequencies  is  very  large.  Some  simplification  approaches  is  intro¬ 
duced,  such  as  the  KMPC  and  A:MPF  approaches.  Results  given  in  Chapter  5 
show  that  the  performance  for  ATMPC  and  fcMPF  is  decent  when  the  spread¬ 
ing  factor  is  high.  When  the  system  load  get  higher,  the  two  approaches  do 
not  perform  well  since  high  MUI  is  introduced.  To  improve  the  performance 
of  the  simplified  metric  in  the  first  iteration,  NSMP  approach  is  introduced 
in  Chapter  4.  Results  display  that  the  performance  of  NSMP  in  the  first 
Iteration  approaches  the  performance  of  the  optimal  metric.  The  reduced 
MUI  in  the  first  iteration  helps  detection  in  later  iterations.  Thus,  the  HER 
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converges  to  a  lower  level  when  NSMP  is  applied  in  the  first  iteration. 

The  most  important  contribution  in  Chapter  5  is  that  the  reason  for  the  error 
floor  is  given.  Basically,  two  factors  lead  to  the  error  floor.  First,  the  ran¬ 
domly  generated  RSC’s  are  not  guaranteed  to  be  uniquely  decodable.  When 
spreading  factor  is  high,  the  problem  is  not  serious  due  to  the  low  proba¬ 
bility  that  two  paths  through  the  super-trellis  have  the  same  code  sequence. 
In  such  circumstance,  the  randomly  generated  codes  perform  well.  When 
spreading  factor  is  close  to  1  or  lower  than  1,  the  randomly  generated  codes 
will  result  in  high  error  floor,  since  the  codes  in  such  condition  are  always  not 
uniquely  decodable.  Note  that  this  part  of  error  floor  exists  for  both  MLSD 
and  iterative  detection. 

The  second  factor  for  the  error  floor  is  the  suboptimality  of  the  iterative 
detection.  In  the  iterative  detection,  though  the  users  borrow  the  a  priori 
information  from  each  other  during  the  iteration,  there  is  no  point  in  the 
consensus  decoder  that  consider  aU  the  users  as  a  whole.  In  this  sense, 
the  consensus  decoder  is  just  partly  cooperative  among  the  users.  (The 
MLSD  is  a  full  cooperative  approach.)  Moreover,  the  flrst  iteration  will 
lead  to  poor  estimation  for  the  symbols,  since  no  apriori  information  is 
available  in  the  flrst  iteration.  When  the  imprecise  a  priori  information  is 
exchange  to  other  users,  MUI  is  introduced  to  other  users.  When  the  users’ 
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codes  are  closely  correlated  with  each  other,  the  MUI  becomes  much  more 
serious.  If  /irMPC  or  A:MPF  is  used  for  simplification,  performance  in  the 
first  iteration  is  deteriorated  further.  Heavy  loaded  system  condition  will  also 
worsen  the  performance  in  the  first  iteration.  The  high  MUI  resulted  from 
the  first  iteration  will  keep  on  affecting  the  performance  of  later  iterations. 
The  reduction  of  the  MUI  in  each  iteration  depends  on  the  error  correction 
capability  of  the  channel  codes.  If  the  channel  codes  are  not  strong  enough 
in  resolving  the  MUI,  high  level  error  floor  is  resulted. 

Prom  results  m  Chapter  5,  we  can  conclude  that  to  lower  the  error  floor 
IS  effective  in  improving  the  system  performance.  To  address  the  problem 
introduced  by  the  non-uniquely  decodable  codes  when  spreading  factor  is  low, 
the  UDBC  codes  axe  introduce  in  Chapter  6.  The  results  show  that  good 
performance  can  be  achieved  by  using  Chang  and  Wolf’s  non-binary  block 
codes  given  in  [17].  However,  we  find  that  the  error  floor  is  still  extremely 
high  if  too  many  users  access  to  the  system.  Thus,  high  coding  rates  cannot 
be  achieved,  which  are  approaching  the  capacity  bound.  The  reason  behind 
this  is  the  suboptimality  of  iterative  detection,  since  we  find  that  the  rate 
1/2  RSC  is  not  strong  enough  to  reduce  the  MUI. 

To  solve  problem  caused  by  the  the  supoptimality  of  iterative  detection, 
several  approaches  can  be  taken.  At  first,  if  PfMPC  or  kUPF  is  used  for 
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simplifying  the  metric,  then  NSMP  should  be  used  in  the  first  iteration  for 
reducing  the  MUI.  Second,  the  channel  interleaver  can  be  used  to  reduce  the 
correlation  among  the  users’  codes.  Thus,  the  MUI  caused  in  each  iteration 
can  be  reduced.  Results  in  Chapter  5  have  shown  the  improvement  by  using 
interleavers.  Since  interleavers  also  increase  the  memory  of  the  super-trellis 
code  of  the  whole  system,  the  performance  given  by  interleavers  is  even  better 
than  the  performance  of  MLSD  without  interleaver,  though  such  comparison 
is  not  fair.  Since  the  application  of  channel  interleavers  does  not  increase  the 
complexity  of  the  decoder  structure,  it  is  an  effective  approach  for  improve 
the  system  performance. 

By  comprehensively  using  the  approaches  discussed  above,  we  achieved  sev¬ 
eral  high  coding  rates  at  the  six  points  in  Figure  7.1.  Performance  at  point  5 
and  6  are  given  in  Chapter  6  by  a  synchronous  system.  At  point  1,  2  and  3, 
both  synchronous  and  asynchronous  system  can  achieve  decent  performance. 
Good  performance  at  point  4  can  only  be  given  by  synchronous  system.  The 
details  of  these  results  can  be  found  in  Chapter  5.  Note  that  all  the  six  points 
are  around  or  higher  than  the  capacity  of  single-user  CDMA.  For  the  first 
four  points,  the  rates  are  equal  to  or  higher  than  the  capacity  of  multiuser 
CDMA.  This  implies  us  that  CFSK  system  has  higher  potential  than  CDMA 
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Capactt/  (bits  pdf  chip) 


CHAPTER  7.  CONCLUSIONS  AND  FUTURE  WORK 


132 


10’ 


Capacity  of  CDMA  and  C-FSK  as  a  function  of  ttia  spreacfing  faca 

C-FSK  with  Intensity  infamaHoti 
Muniuser  CDMA 

—  C~FSK  without  intensity  Information 
Conventional  CDMA 


10" 


lo'T 


10’ 


10 


1.  Rate  1/6  code,  64ftequencies,12us£s 

2.  Rate  2/6  code,64  &equenctes,12  useis 

3.  Rate  1  /3  code,8  fx)equetic!tes,4t;)>se]:s 

4.  Rate  1  /2  code,  4  fioequencdes,  3  users 

5.  Rate  1  /2  code,  3  ftequencies,  6  users 

6 .  Rate  1  /2  co de,3  frequencie s,l 4  users 


10’ 


Spreading  Factor  N/K  (chips  per  irser) 


Figure  7.1:  Achieved  sum  rates  compared  with  capacity  bounds. 
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7.2  Future  Work 

A  problem  caused  by  the  signal  mapping  in  our  scheme  is  that  it  reduces  the 
free  distance  of  the  convolutional  codes.  Performance  of  the  convolutional 
codes  is  thus  degraded  and  the  codes  are  not  strong  enough  to  reduce  the 
MUI  when  a  large  number  of  users  access  to  the  system.  Possible  approach  to 
improve  the  performance  of  the  users’  component  decoders  is  to  map  the  bit 
sequence  into  the  frequency  sequence  by  non-binary  codes  with  decent  dis¬ 
tance  performance.  Nonbinary  codes,  such  as  convolutional  codes  over  GF{q) 
[41]  [41]  [42]  [40]  or  convolution  codes  over  rings  [43]  [44]  [45]  [46]  [47],  are  candi¬ 
dates  for  substituting  the  binary  convolutional  codes  in  the  former  system. 
Using  these  nonbinary  codes  as  building  blocks,  the  more  powerful  nonbi¬ 
nary  Turbo  codes  can  be  constructed  [48]  [49].  The  corresponding  decoding 
algorithms  [50]  [51]  can  be  adapted  for  applications  in  CFSK  system. 

For  the  randomly  generated  nonbinary  codes,  unique  decodability  of  these 
codes  is  still  a  major  problem  for  applications  of  low  spreading  factors.  To 
solve  the  problem,  uniquely  decodable  block  codes  with  large  codeword  num¬ 
bers  are  useful.  The  KS  code  given  by  [37]  and  [39]  provides  us  some  inspi¬ 
ration  for  using  the  Reed-Solomon  code  to  construct  the  imiquely  decod¬ 
able  block  code  set.  Also,  researches  of  other  nonbinary  superimposed  codes 
[37]  [38]  [39]  also  provide  us  clues  for  constructing  the  uniquely  decodable  block 
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code  set. 

We  may  also  avoid  generating  the  codes  randomly.  Instead,  constructions  for 
the  nonbmary  uniquely  decodable  trellis  codes  can  be  used.  In  this  approach, 
we  need  to  extend  former  results  for  the  uniquely  decodable  trellis  codes  on 
the  binary  adder  channel  [26]  [52]  to  nonbinary  uniquely  decodable  trellis 
codes  on  the  CPSK  channel  with  intensity  information. 

For  future  application  of  the  CFSK  system,  performance  of  CFSK  system  on 
fading  channels  and  the  resistance  of  CFSK  system  to  the  near-fax  effects  need 
to  be  studied  in  detail.  Though  the  CFSK  system  is  supposed  to  have  decent 
performance  under  these  channel  conditions,  a  theoretical  investigation  is 

indeed  necessary  to  prove  the  effectiveness  of  the  CFSK  system  in  dealing 
with  these  problems. 
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